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FOREWORD 


This  report  is  one  of  two  computer  program  user's  manuals  prepared 
by  Aerotherm  Corporation  under  USAF  Contract  F04611-70-C-0012.  Included  herein 
is  Volume  I  of  the  manual  for  Version  3  of  the  Aerotherm  Charring  Material 
Ablation  code.  This  volume  describes  the  problems  solved  by  the  code  and  pre¬ 
sents  an  input  (card  format)  user's  guide  and  sample  problems.  The  report  was 
first  published  as  Aerotherm  Report  No.  UM-70-14.  The  work  was  administered 
under  the  direction  of  the  Air  Force  Rocket  Propulsion  Laboratory,  Motor  Compo¬ 
nent  Development  Branch  with  Mr.  R.  J.  Schoner  as  project  officer. 

Mr.  M.  R.  Wool  was  program  manager  and  principal  investigator. 
Significant  additional  assistance  was  also  provided  by  Dr.  C.  B.  Moyer. 

This  technical  report  has  been  reviewed  and  is  approved. 


R.  J.  Schoner 
Project  Engineer,  AFRPL 


ABSTRACT 


This  two-volume  report  describes  a  Fortran  IV  computer  code  which 
computes  the  transient  thermal  and  ablation  response  of  a  charring  ins -ration 
material  structure.  The  program  is  for  one-dimensional  bodies,  but  can  treat 
a  variety  of  shapes,  including  planes,  cylinders,  spheres,  and  more  general  • 
thermal  “stream  tube"  bodies.  The  program  can  treat  complex  systems  including 
a  main  ablating  material,  several  charring  back-up  materials,  and  a  multiple 
non-charring  material  back-up  structure. 

An  unusual  feature  of  the  code  is  the  very  general  heated  surface 
boundary  condition,  which  can  account  for 

•  Simple  specified  temperature  and  recession  rate 

•  Specified  heat  flux  with  no  recession 

•  General  thermochemical  erosion  model  incorporating  complete 
chemical  erosion  computations,  both  equilibrium  and  non- 
equilibrium,  for  any  material  exposed  to  any  environment 

The  code  has  seen  extensive  use  for  thermal  performance  studies  of 
ablating  heat  shields,  rocket  nozzles,  and  spacecraft  structures. 

Volume  I  of  this  report  contains  descriptions  of  the  problem  treated 
the  equations  solved,  the  input  information  required  of  the  program  user,  and 
the  program  output  information.  It  also  provides  a  card-by-card  user's  input 
guide  and  a  number  of  sample  problem  input  and  output  listings.  Volume  II  of 
the  report  contains  supplemental  information  on  the  specific  Fortran  IV  codings 
It  includes  program  listings,  flow  charts,  and  definitions  of  Fortran  variable 


names. 


When  U.S.  Government  drawings,  specifications,  or  other  data 
are  used  for  any  purpose  other  than  a  definitely  related 
Government  procurement  operation,  the  Government  thereby  Incurs 
no  responsibility  nor  any  obligation  whatsoever,  and  the  fact 
that  the  Government  may  have  formulated,  furnished,  or  in  any 
way  supplied  the  said  drawings,  specifications,  or  other  data, 
is  not  to  be  regarded  by  implication  or  otherwise,  or  in  any 
manner  licensing  the  holder  or  any  other  person  or  corporation 
or  conveying  any  rights  or  permission  to  manufacture,  use,  or 
sell  any  patented  invention  that  may  in  any  way  be  related 
thereto. 
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SECTION  1 


INTRODUCTION  AND  GENERAL  PROGRAM  DESCRIPTION  / 

1.1  PURPOSE  OF  THIS  DOCUMENT 

This  report  is  Volume  I  of  a  two  volume  user's  description  of  the 
Aerotherm  Charring  Material  Thermal  Response  and  Ablation  computer  program. 
Version  3  (CMA) .  It  has  as  its  chief  goals  the  description  of  the  preparation 
of  the  necessary  input  data  for  the  program  and  the  specification  of  various 
operational  details  of  the  program  such  as  array  sizes  and  tape  numbers  used. 

As  supplementary  information,  this  volume  includes  a  very  brief  description 
of  the  computations  performed.  This  description  is  not  intended  to  be  com¬ 
plete,  but  merely  serves  to  give  the  user  some  feel  for  the  calculations  as 
well  as  a  better  understanding  of  the  input  requirements. 

Users  interested  in  a  more  complete  description  of  the  program  itself 
are  referred  to  Reference  1.  An  earlier  version  of  this  program  is  described  in 
detail  in  Reference  2.  This  user's  manual  pertains  specifically  to  what  will 
be  termed  Version  3  of  the  CMA  code.  It  differs  from  the  earlier  Version  2 
edition  described  in  Reference  3  chiefly  in  the  added  capability  to  treat 
charring  back-up  materials.  Hence,  an  alternate,, designation  '  (CMA/CBM  or  CBM)— 
is  used  occasionally  in  this  manual.  '' 

Volume  II  of  this  user's  manual  contains  detailed  information  about 

► 

Fortran  IV  codings.  These  include  Definitions  of  Fortran  variables.  Flow  Charts 
of  Fortran  routines,  and  Listings  of  Fortran  routines. 


1.2  GENERAL  DESCRIPTION  OF  THE  CMA  PROGRAM 

The  Charring  Material  Thermal  Response  and  Ablation  Program,  Version  3, 
is  an  implicit,  finite-difference  computational  procedure  for  computing  the  one¬ 
dimensional  transient  transport  of  thermal  energy  in  a  three-dimensional  iso¬ 
tropic  material  which  can  ablate  from  a  front  surface  and  which  can  decompose 
in-depth.  Decomposition  reactions  are. based  on  a  three  component  model.  The 
program  permits  up  to  eight  different  back-up  materials  of  arbitrary  thickness, 
five  of  which  may  char.  The  back  wall  of  the  composite  material  may  transfer 
energy  by  convection  and  radiation. 


I 


The  ablating-surface  boundary  condition  may  take  one  of  three  forms: 

Option  1  -  General  convection-radiation  heating  with  coupled 

mass  transfer,  using  a  transfer  coefficient  approach, 
including  the  effects  of  unequal  heat  and  mass  trans¬ 
fer  coefficients  (non-unity  Lewis  number)  and  unequal 
mass  diffusion  coefficients.  Surface  thermochemistry 
computations  need  not  presume  chemical  equilibrium  at 
the  surface,  hut  can  account  for  kinetic  effects,  and 
can  also  consider  thin  layer  melting  or  failing. 

Option  2  -  Specified  surface  ^temperature  and  surface  recession 

rate.  / 

Option  3  -  Specified  radiation  view  factor  and  incident  radia¬ 

tion  flux,  as  functions  of  time,  for  a  stationary 
surface.  ' 

Any  combination  of  options  may  be  used  for  a  single  computation.  Op¬ 
tion  3  is  appropriate  to  cooldown  after  termination  of  convective  heat  input  and 
is  often  useful  in  conjunction  witu  Options  1  and  2. 

The  following  sketch  illustrates  the  general  physical  problem  treated 
by  the  CMA  code.  - 
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As  the  material  is  heated,  the  original  virgin  material  (or  rather  one  or  more 
components  of  the  original  composite  virgin  material)  pyrolyzes  and  yields  a 
pyrolysis  gas,  which  percolates  away  from  the  pyrolysis  zone,  and  a  porous  resi¬ 
due,  which  for  most  materials  of  interest  is  a  carbonaceous  char,  possibly  re¬ 
inforced  with  refractory  fibers  or  cloth. 


Superimposed  on  this  basic  problem  may  be  a  number  of  even  more  com¬ 
plex  events.  The  pyrolysis  gases  percolating  through  the  char  may  undergo 
further  chemical  reactions  among  tuemselves ,  and  may  react  with  the  char,  either 
eroding  it  or  depositing  additional  residue  upon  it  (“coking").  The  char  itself 
may  collapse  or  fragment  from  mechanical  or  thermal  stresses,  and  the  refractory 
reinforcements  may  melt  or  suffer  mechanical  damage.  Finally,  various  consti¬ 
tuents  of  the  residue  structure  may  react  chemically  with  each  other,  changing 
the  nature  of  the  char,  and  various  mechanical  forces  may  remove  material  from 
the  surface. 

Despite  these  complexities,  it  is  found  that  the  "simple  physics" 
described  by 


virgin  plastic  -*■  char  +  gas 

underlies  a  wide  range  of  problems  of  technical  interest,  and  for  a  great  many 
materials,  such  as  carbon  phenolic,  graphite  phenolic,  and  wood,  constitute  all 
the  events  of  interest.  Such  events  as  coking,  mechanical  erosion,  melting, 
and  subsurface  reactions  (other  than  pyrolysis)  are  less  common  curd  generally 
characterize  specific  problems. 

The  CMA  Code,  Version  3,  treats  the  basic  pyrolysis  or  charring  event 
for  a  main  insulation  material.  This  material  may  be  backed-up  by  a  number  of 
charring  back-up  materials,  which  in  turn  may  be  backed-up  by  noncharring  struc¬ 
tural  elements.  The  CMA  code  will  treat  cracking  or  fissuring  of  the  char,  as 
well  as  surface  mechanical  removal  or  melting  of  the  char  structure. 
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SECTION  2 


PROGRAM  FEATURES,  CAPABILITIES,  AND 
COMPUTATIONAL  PROCEDURES 

2.1  IN-DEPTH  SOLUTION  ASPECTS 

The  in-depth  solution  procedure  is  the  controlling  feature  of  the  code 
and  accounts  for  most  of  the  computing  time.  This  procedure  is  basically  a 
transient  heat  conduction  calculation  which  is  coupled  to  a  pyrolysis  rate  cal¬ 
culation  (and  to  boundary  conditions  discussed  in  Section  2.4  below).  The  fol¬ 
lowing  subsections  discuss  important  aspects  of  the  in-depth  solution. 

2.1.1  Decompos i ton  (Pyrolysis  or  Charring) 

Since  many  decomposing  char  forming  materials  appear  to  behave  as 
three  independently  pyrolyzing  components,  the  program  uses  a  three-component 
decomposition  model  for  the  main  material  and  for  any  decomposing  back-up 
material. 

The  resin  filler  is  presumed  to  consist  of  two  components  which  de¬ 
compose  separately,  while  the  reinforcing  material  is  the  third  component  which 
can  decompose.  The  instantaneous  density  of  the  composite  is  given  by 

p  =  r(oA  +  pb),  +  (i  -  r)pc  (i) 

where  A  and  B  represent  components  of  the  resin,  and  C  represents  the  reinforc¬ 
ing  material,  T  is  the  volume  fraction  of  resin  and  is  an  input  quantity  (Sec¬ 
tion  3.1.2).  Separate  r  values  may  be  entered  for  any  decomposing  back-ups. 

Each  of  the  three  components  can  decompose  following  the  relation 


(2) 


where  pr.  is  the  residual  or  terminal  density  of  component  i,  and  p_.  is  the 
x  ui 

original  density  of  component  i.  The  values  pQ^  pr^,  B^,  and  Ea.  are  input 

parameters  for  i  =  A,B,C  for  the  main  material.  A  separate  i  =  A,B,C  set  of 

pyrolysis  data  may  be  entered  for  each  decomposing  back-up. 
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The  decomposition  event  is  computed  explicitly.  That  is,  "old"  tem¬ 
peratures  are  used  in  Equation  (2)  to  compute  the  decomposition  rates  and  the 
new  nodal  densities. 


i 


I 


2.1.2  In-Depth  Energy  Balance 

For  the  purpose  of  writing  the  in-depth  energy  balance  differential 
equation,  we  introduce  the  x  coordinate  system  tied  to  the  receding  surface,  as 
shown  in  the  following  sketch.  In  this  system,  we  have  the  energy  equation  of 
the  form 


PC 


3T\ 
P  36  )y 


1  3 _ 

A  3x 


+ 


+  5  PCp 


m  3h  \ 

+  -2.  —2.) 
A  3x  JQ 


(3) 


?n  which  the  individual  terms  have  physical  meanings  which  may  be  interpreted 
as  follows  (from  left  to  right):  rate  of  storage  of  sensible  energy,  net  rate  of 
thermal  conduction,  pyrolysis  energy  rate,  convection  rate  of  sensible  energy 
due  to  coordinate  system  movement,  and  net  rate  of  energy  convected  with  pyroly¬ 
sis  gas  passing  a  point. 


In  this  equation,  the  local  specific  heat  and  thermal  conductivity 
are  formulated  from  input  temperature  functions  for  both  virgin  plastic  and 
cher,  C  ,  C  ,  k^ ,  and  k  (see  Section  3.1.7  below).  In  partially  pyrolyzed 

PP  pC  r  C 

sones  (p  <  p.  <  o  ) ,  the  specific  heat  is  formulated  with  a  special  mixing  rule 


:  =  xC  +  (1  -  x)C 

P  PP  pc 


(4) 


where  the  weighting  variable  x  is  based  on  the  convenient  fiction  that  partially 
pyrolyzed  material  is  a  simple  mixture  of  pure  virgin  plastic  and  pure  char. 
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The  quantity  x  is  defined  as  the 
mixture  which  yields  the  correct 

x  = 

The  thermal  conductivity  k  will  be  weighted  in  the  same  manner  unless 
the  user  provides  specific  information  on  a  more  complex  x  weighting  in  the 
input  data.  That  is,  the  user  may  input  functions  of  x  for  the  following  equa¬ 
tion  for  the  thermal  conductivity  of  partially  degraded  plastic: 

k  -  f1M*p  +  f2(x)lcc  (6) 

References  4  and  5  present  an  account  of  the  rationale  for  the  use 
of  f-functions,  along  with  recommended  values. 

ZT 

The  pyrolysis  gas  enthalpy  is  an  input  temperature  dependent  function 
(see  Section  3.1.8  below).  The  quantity  h  =  (pphp  “  hcpc^^pp  ~  Pc)  "''s  coinPuted 
from  temperature  dependent  hp  and  hc  values  determined  from  the  input  tempera¬ 
ture  dependent  Cpp  and  Cpc  values  and  input  enthalpies  of  formation  AHfp  and 

AHfc* 

Determination  of  the  local  cross-section  area  A  will  be  discussed  in 
Section  2.1.3  below. 

The  other  quantities  in  Equation  (3) ,  T,  p  and  S  are  dependent  quan¬ 
tities  discovered  during  the  solution  process  (except  that  in  some  boundary 
condition  options  S  may  be  input) . 

2.1.3  Nodal  Coordinate  Layout  and  Body  Shape  Specification 

The  basic  solution  procedure  of  the  CMA  code  is  of  the  finite  differ¬ 
ence  type.  The  following  sketch  illustrates  a  typical  material  section  that  is 
to  be  analyzed.  Nodal*  positions  are  specified  by  tabulating  the  thickness  of 
successive  nodes  (see  Section  3.1.5  below).  The  program  permits  a  thermal  con¬ 
tact  resistance  between  adjacent  nodes,  and  thus  between  adjacent  materials 
(see  Section  3.1.5  below).  Nodes  are  numbered  sequentially  from  the  surface. 

A  total  of  100  nodes  is  permitted. 


mass  fraction  of  pure  plastic  in  this  imaginary 
local  density: 


* 

As  used  here,  the  term  node  refers  to  a  finite  region  in  the  slab  wall,  not  to 
a  discrete  point  or  plane. 
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Ir.  harmony  with  the  shifting  x-coordinate  system  introduced  in  Sec- 
tion  2.1.2,  nodal  coordinates  are  tied  to  the  heated  surface.  Any  surface  reces¬ 
sion  causes  the  last  (deepest)  node  in  the  ablating  material  to  shrink.  If  sur¬ 
face  recession  proceeds  far  enough,  this  node  is  dropped  from  the  rear  surface 
of  the  ablating  material  and  the  number  of  nodes  reduced  by  one.  Computation 
ceases  when  only  one  node  of  the  main  material  remains. 


Ablating  surface 


•  n  =  2 


Last  node  in  main  ablation  material 
(consumed  as  surface  recedes) . 

First  node  in  backup  material 


_ • _  Last  node  in  network 

Back  wall 

Sketch  of  Nodal  Network 

The  nodes  are  used  for  all  energy  balance  calculations. 

Since  decomposition  rates  depend  exponentially  on  temperature  the 
decomposition  calculations  in  the  main  material  proceed  on  the  basis  of  a 
"nodelet"  network  using  small  nodes  a  fraction  as  large*  as  the  specified  nodal 
layout  used  in  energy  balance  calculations.** 

The  decomposition  events  are  computed  for  each  nodelet.  The  total 
density  change  for  each  node,  is  taken  as  the  sum  of  the  nodelet  density  change 
rates . 

The  program  permits  the  specification  of  a  number  of  geometries: 

1.  Plane 

/* 

2.  Cylindrical  or  annular,  with  heated  surface  either  inner  or  outer 

Usually  one-tenth,  although  the  number  of  nodelets  per  node  may  be  selected 
by  the  user  (see  Section  3.1.3  below). 

**The  nodelet  device  is  not  used  in  charring  back-up  materials,  which  generally 
pyrolyze  rather  slowly. 


3.  Spherical  or  spherical  shell,  with  heated  surface  either  inner 
or  outer 

4.  General  "thermal  stream  tube"  geometry,  area  varying  as  depth 
to  any  power 

5.  General  "thermal  stream  tube"  geometry,  area  varying  arbitrarily 
with  depth. 

The  distinction  between  these  five  options  is  indicated  by  the  user  according 
to  the  input  rules  given  in  Section  3.1.5  below. 

2.1.4  Computation  Procedure 

2 . 1 . 4 . 1  General  Pattern 

As  in  all  finite  difference  procedures,  each  step  of  the  solution  is 
made  over  an  incremental  time  step  A0 .  Each  computational  step  has  three  main 
events:  internal  decomposition,  internal  nergy  balance,  and  surface  boundary 
energy  balance .  Computation  c-r  the  thrr  i  events  gives  "new"  values  of  nodal 
densities,  pyrolysis  gas  flow,  nodal  temperatures,  and  surface  temperature  and 
ablation  rate.  The  program  is  then  ready  for  the  next  step. 

The  pyrolysis  event  is  computed  for  each  node let  in  the  main  material 
from  Equation  (2),  summed  over  the  three  components.  Summation  over  all  nodelets 
in  a  node  gives  the  nodal  3 p/3 8  and  contribution  to  mg.  Note  that  "old"  known 
nodal  temperatures  (interpolated  into  the  nodelet  locations)  are  used  in  the 
pyrolysis  calculation. 

The  internal  energy  balance  equation  then  is  computed  "implicitly" 
for  each  node,  using  "new"  temperatures  in  the  heat  conduction  terms.  The  en¬ 
ergy  balance  is  linked  explicitly  to  the  decomposition  events,  however,  since 
the  pyrolysis  gas  fluxes  used  in  the  energy  balance  are  derived  from  the  explicit 
decomposition  calculation.  The  energy  balance  is  also  linked  explicitly  to  the 
surface  bound lry  condition  through  the  use  of  an  "old"  surface  recession  rate  in 
all  convection  terms  involving  fluxes  of  solids.  All  other  links  to  surface 
events  are  imp licit. 

The  "implicit"  formulation  of  the  in-depth  nodal  energy  equations 
yields  a  tri-diagonal  set  of  equations  which  is  solved  for  the  unknown  tempera¬ 
tures  in  two  passes  of  direct  elimination.  The  first  pass  eliminates  one  un¬ 
known  from  each  equation  and  leaves  the  equation  for  the  first  node  (at  the 
heated  surface)  with  only  one  unknown.  In  surface  boundary  condition  Option  2, 
this  temperature  is  specified  and  so  is  known.  In  Options  1  and  3,  the  first 
node  equation  is  non-linear  and  so  the  first  nodal  temperature  must  be  discovered 
by  an  iterative  procedure.  This  will  be  discussed  in  Sections  2.4.2  and  2.4.3. 
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In  any  case,  once  this  nodal  temperature  is  determined,  the  second  elimination 
pass  determines  the  other  unknown  nodal  temperatures  and  the  entire  in-depth 
calculation  is  ready  for  a  new  time  step. 

2. 1.4. 2  Limitations  on  Time  Step  Size 

Experience  shows  that  the  implicit  solution  procedure  of  the  CMA  pro¬ 
gram  is  almost  always  stable  regardless  of  time  step  size,  and  thus  has  an  im¬ 
portant  advantage  over  explicit  procedures  for  which  stability  considerations 
severely  limit  the  allowable  time  step  size. 

However,  in  order  to  assure  a  smooth  progression  of  the  solution, 
limitations  must  be  imposed  upon  the  time  increments.  The  user  may  specify 
as  input  (see  Section  3.1.3  below)  a  maximum  allowable  time  step  not  to  be 
exceeded  under  any  circumstances .  Certain  other  limits  are  computed  automati¬ 
cally  and  may  reduce  the  solution  time  step  A0  below  this  value: 

•  Time  steps  are  sized  to  limit  the  change  of  surface  temperature 
during  one  step  to  approximately  50°R 

•  At  the  initiation  of  a  problem,  or  when  discontinuities  are  noted 
in  the  surface  boundary  condition  option,  the  time  step  is  re¬ 
duced  to  0.01  seconds,  with  subsequent  values  limited  to  twice 
the  preceding  value  of  the  time  increment. 

•  Time  steps  are  limited  so  that  recession  during  a  time  step, 

SA0 ,  does  not  exceed  the  thickness  of  the  smallest  nodelet  in 
the  system. 

•  After  all  these  limits  have  been  applied,  the  time  step  is  re¬ 
duced  so  that  the  time  remaining  to  the  next  printout  txme  is 
divided  into  an  integral  number  oir  time  steps  of  the  current 
size  (so  that  no  sudden  reduction  in  time  step  is  required  in 
order  to  match  a  printout  time) . 

2.2  INITIAL  CONDITION 

As  initial  conditions,  the  user  may  specify  the  temperature  of  each 
•ode  and  in  addition  may  call  out  each  node  in  a  charring  material  as  either 
virgin  plastic  or  pure  char.  This  allows  studies  of  pre-charred  materials  Ana 
re-start  problems.  Section  3.1.5  below  describes  the  necessary  input  rules. 

2.3  BACKWALL  BOUNDARY  CONDITION 

Heat  transfer  from  the  backwall  of  the  system  is  treated  with  a  sim¬ 
ple  temperature-potential  convective  transfer  plus  a  radiative  term,  both 


1 


,  .  > 

communicating  with  a  “reservoir1*  at  temperature  ^Ies-  The  convective  coeffi¬ 
cient,  backwall  emissivity  (or  emittance)*,  and  Tres  are  all  input  constants 
(see  Section  3.1.5  below). 

t. 

2.4  ABLATING  SURFACE  BOUNDARY  CONDITIONS 

The  ablating  surface  boundary  conditions  may  take  one  of  three  forms 
at  any  instant,  as  determined  by  the  user  in  the  make-up  of  the  time  dependent 
boundary  condition  input  table  (see  Section  3.1.9  below).  The  two  simpler  op¬ 
tions  will  be  described  first,  even  though  for  historical  reasons  these  are  de¬ 
noted  "Option  2"  and  "Option  3".  The  general  thermochemical  ablation  boundary 
condition,  denote  1  "Option  1",  is  discussed  last. 

2.4.1  Option  2,  Specified  Surface  Temperature  and  Surface  Recession  Rate 

/ 

In  this  option,  surface  temperature  and  recession  rat^  are  specified 
by  the  User  as  input  f mictions  of  time.  All  surface  energy  balance  and  surface 
thenoochemical  considerations  are  bypassed.  Section  2.1. 4.1  above  describes 
how  this  specified  temperature  can  be  inserted  into  the  in-depth  solution  pro¬ 
cedure  after  the  first  pass  of  direct  elimination. 

2.4.2  Option  3,  Specified  Input  Heat  Flux  with  Zero  Surface  Recession 
Rate  (Cooldown  Option) 

The  following  sketch  illustrates  the  surface  energy  balance  situation 
for  this  option.  Equation  (7)  is  the  appropriate  surface  energy  balance  equa¬ 
tion 

^rad  ^jrad 

Jin  jout 

^cond 


'  a  q  —  F_cjc  T*  -  q  =  0  (7) 

w^rad  3  w  w  scond  '  7 

In  this  equation,  aw  is  the  surface  absorptance  and  is  assumed  to  equal  ew, 
the  surface  emissivity,  which  is  computed  from  input  temperature  dependent  prop¬ 
erty  tables  (see  Section  3.1.7  below).  If  the  surface  material  is  partially 
pyrolyzed,  ew  is  computed  according  to  the  mixture  rule  of  Equation  (4)  .  The 
term  qra(j  is  an. input  time  dependent  heat  flux  specified  in  the  boundary  condi¬ 
tion  time  tables  (see  Section  3.1.9  below) .  F^  is  an  input  view  factor  for 
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Option  3  calculations.  It  may  depend  on  time  (see  Section  3.1.9  below).  The 
term  <lcon(j  represents  tne  rate  of  heat  conduction  into  the  material.  It  is  the 
"mechanism"  which  joins  the  surface  energy  balance  equation  to  the  first  nodal 
energy  balance  (since  qconcl  appears  in  the  nodal  energy  balance  equation)  and 
thus  to  the  in-depth  solution  procedure,  as  discussed  in  Section  2. 1.4.1  aboVe. 
The  first  pass  of  the  elimination  process  fot  the  in-depth  tri-diagonal  system 
reduction  leaves  the  nodal  energy  balance  equation  for  the  first  node  as 


%on«  ‘  +  B  !») 

Equation  (8),  when  substituted  into  the  surface  energy  balance  equation  (7), 
yields  a  non-linear  equation  in  Tw  which  the  CMA  program  solves  iteratively  by 
the  Newton- Raphson  method. 


2.4.3 


Option  1,  General  Convective  Heating  and  Thermochemical 
Erosion  Option 


2. 4. 3.1  General  Description  and  Basic  Energy  Equation 

In  this  option,  events  at  the  heated  surface  are  determined  by  con¬ 
vective  heating  and  by  surface  thermochemical  interactions  with  the  boundary 
layer  gases.  The  sketch  below  illustrates  the  surface  control  volume  and  the 
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Sketch  of  Surface  Energy  Balance  Control  Volume 
and  Energy  Fluxes  in  Option  1 


energy  fluxes  of  interest.  The  surface  energe  balance  equation  employed  is 
of  the  convective  transfer  coefficient  type.  In  the  CMA/CBM  program,  this 
energy  balance  equation  takes  the  following  form: 


peueCH(Hr  '  hew>  +  pe 


w[£|zie  -  zs»>hi  -  B'h»] +  ichc +  v-c 

-  y 

^chem 


+  “w^rad  "  FoewTw  .  "  %ond  =  0 
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The  last  three  terms  in  Equation  (9)  are  the  same  as  in  the  Option  3 
boundary  condition  previously  discussed  in  Section  2.4.2  above.  The  other 
terms  in  this  equation  require  some  further  comments . 

Before  commencing  a  term  by  term  discussion  of  Equation  (9) ,  however, 
it  will  be  useful  to  describe  the  general  nature  of  this  transfer  coefficient 
expression.  Like  all  such  expressions.  Equation  (9)  is  an  approximation,  the 
usefulness  of  which  depends  mainly  on  the  validity  of  the  transfer  coefficient 
approach.  A  discussion  of  this  subject  As  far  beyond  the  scope  of  the  present 
document.  It  may  be  observed  here  that  transfer  coefficients  have  successfully 
correlated  both  data  and  "exact"  solutions  in  simple  heat  or  mass  transfer  prob¬ 
lems,  and  in  combined  heat  and  mass  transfer  problems  for  unity  (or  near  unity) 
Lewis  number.  Equation  (9)  attempts  to  extend  the  transfer  coefficient  approach 
to  both  non-unity  Lewis  number  and  unequal  mass  diffusion  coefficient  problems, 
still  allowing  for  chemical  reactions  and  net  mass  transfer,,  effects .  This 
approach  was  suggested  in  Reference  (6) .  Its  validity  is  discussed  in  Refer¬ 
ences  1  and  7. 

/ 

In  Equation  (9) ,  the  term  <3sen  represents  the  "sensible  convective  / 
heat  flux."  Physically,  this  is  the  convective  heat  flux  which  would  occur  // 
for  a  frozen  boundary  layer  and  a  non-catalytic  -wall  in  the  absence  of  mass  / 
transfer;*  it  excludes  all  chemical  energy  contributions.  (The  term  q 

sen 

is  perhaps  more  usually  written  in  the  form 


*sen 


='  peUeCH(V 


-  h  ) 


(10) 


but,  since  generally  it  is  more  convenient  for  the  user  to  input  Hr  rather  than 
V  <Jsen  in  Equation  (9)  has  been  written  in  a  modified  form  in  which  Hr  appears. 
This  form  has  tSie  additional  advantage  that  the  driving  force  for  energy  trans¬ 
fer  involves  only  edge  gas  states.  The  derivation  of  the  modified  form  from 
Equation  (1)  is  given  in  Reference  1  and  Reference  6. 

The  transfer  coefficient  PeueCH  and  the  recovery  enthalpy  Hr  are 
time  dependent  variables  input  by  the  ptogram  user  in  appropriate  time  tables 
(see  Section  3.1.9  below).  The  transfer  coefficient  is  automatically  modified 
from  the  input  value  to  account  for  various  effects,  as  discussed  in  Section 
2. 4. 3. 4  below.  The  quantity  h^  is  part  of  the  input  thermochemical [data  dis¬ 
cussed  below. 


n - 

More  generally  in  the  presence  of  chemical  reaction  it  is  the  diffusive  heat 
flux  from  the  gas  to  the  wall  even  in  the  presence  of  net  mass  transfer,  pro¬ 
vided  the  boundary  layer  is  frozen  and  the  wall  is  catalytic. 
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The  term  qnh^m  represents  the  net  of  a  number  of  fluxes  of  chemical 
energies  at  the  surface.  The  Z*-dif ference  term  represents  transport  of  chemi¬ 
cal  energy  associated  with  chemical  reactions  at  the  wall  and  in  the  boundary 
layer;  it  is  the  chemical  energy  parallel  to  the  sensible  convective  heat  flux 
term.  The  2*  driving  forces  for  diffusive  mass  transfer  include  the  effects 
of  unequal  diffusion  coefficients;  for  equal  diffusion  coefficients  the  Z*'s 

reduce  to  ‘■he  familiar  nu_.  ■?  fractions  K.  The  m  h  and  m  h  terms  represent 

c  c  9  9 

energy  fluxes  arriving  at  the  surface  from  within  the  solid  material  and  the 
B'h^  term  represents  energy  leaving  the  surface  in  the  gross  motion  (blowing) 
of  the  gas  adjacent  to  the  surface. 

The  source  of  the  various  terms  in  qc^em  can  better  be  visualized 

after  a  slight  regrouping  of  terms.  Observing  that  B'  =  +  B^,  and  bringing 

the  m  h  term  inside  the  brackets,  we  have 
c  c 
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Of  the  quantities  in  this  qcJiem  expression,  the  convective  mass  transfer 

coefficient  p  u  C  is  obtained  from  the  input  time  dependent  values  of 
c  e  m 

p  ueCH  (after  any  of  the  adjustments  described  in  Section  2.4. 3.4  below),  by 
multiplying  peucCH  by  an  input  constant  value  CM/CH  (see  Section  3.1.10.2  below). 
The  quantity  m  will  have  been  determined  previous  to  the  energy  balance  calcu¬ 
lations  in  each  time  step  in  tiife^decoraposition  section  of  the  code  (see  Section 
,  2. 1.4.1  above) .  The  quantity  B^  may  be  found  as  its  definition  B^  =  mg/PeuecM* 

The  enthalpy  hc  is  obtained  from  integrations  of  the  input  temperature  dependent 
specific  values  Cp  (see  Section  3.1.7  below).  The  pyrolysis  gas  enthalpy  is 
obtained  from  the  input  temperature  dependent  hg  table  (see  Section  3.1.8  below). 

Remaining  quantities  not  yet  discussed  are  B^,  Tw  (which  does  not 
appear  explicitly  but  which  is  necessary  to  evaluate  the  temperature  dependent 

values  of  the  quantities  hc  and  hg) ,  h^,  X)z£ehiW»  £Z£whiW'  and  hw'  The 

quantities  T  ,  £z*whiW'  and  hw  are  inPut  by-  the  user  as  the  dependent  variables 

in  a  table  with  three  independent  variables:  P,  B',  and  B'.  Similarly,  the 
T  “  c 

quantities  ]Cz*ehiW  and  he  are  inPut  as  the  dependent  variables  in  a  table 

with  P  and  T  as  independent  variables.  These  tables  are  typically  generated 

by  the  thermochemistry  codes  described  in  Section  2. 4. 3. 2  below.  Further 

discussion  of  these  tables  is  given  in  Section  3.1.10.3. 
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The  code  surface  energy  balance  solution  procedure  may  be  summarized 

as  follows: 

1.  Look  up  Hr,  PeueCjj/  P,  and  qra(J  in  input  functions  of  time  tables 

2.  Correct  or  adjust  PeueCjj  for  various  effects 

3.  Compute  PeueCM  =  C^/C^ (peueCH) 

4.  Compute  B*  =  m  /p  u  C„ 

g  g'  *^e  e  M 

5.  Assume  B' 

c 

6.  With  P,  Bg,  B£,  look  up  in  input  surface  thermochemistry  tables 
values  of  Tw,  Ez*whiW'  hw 

7.  With  p  and  ^  ^ook  up  in  input  edge  gas  thermochemistry  table 
values  of  £s{ehiwV  he 

8.  With  Tw,  look  up  in  input  property  tables  values  of  ew,  hc,  h^ 

9.  Obtain  values  of  A  and  B  in  expression  cjcon(j  =  ATw  +  B  from  in- 
depth  nodal  energy  balance  solution  routine 

10.  Construct  Equation  (9),  noting  departure  from  zero,  if  any 

11.  Adjust  B^  guess  to  reduce  departure  from  zero 

12.  Go  to  Step  (6)  and  continue 

This  procedure  converges  on  a  new  B^.  value  in  very  few  iterations.  The  same 
procedure  is  used  with  Tw  as  the  independent  variable  and  as  a  dependent 
variable  in  portions  of  the  table  where  B^  varies  slowly  or  not  at  all  with  tw 
(see  Section  3.1.10.3  below  for  input  rules  and  suggestions). 

2. 4. 3. 2  Use  of  Thermochemistry  Codes  to  Generate  Input  Data 

Section  2. 4. 3.1  above  makes  it  clear  that  the  CMA  code  requires  some 
complex  tabular  thermochemical  input  if  the  Option  1  boundary  condition  is  to  be 
used.  These  tables  are  generated  by  any  one  of  a  number  of  separate  computer 
codes.  The  most  recent  such  code  is  designated  the  Equilibrium  Surface  Thermo¬ 
chemistry  Code,  Version  3  (EST3) .  It  is  a  general  open  and  closed  system  ther- 
mochemical  equilibrium  code  specifically  constructed  for  this  purpose.  It  is 
described  in  Reference  8.  A  generally  similar  code  which  differs  from  EST3 
only  in  added  detail  is  designated  ACE  and  is  described  in  Reference  9.  An  • 
older  version  of  EST3  was  designated  EST2  and  is  described  in  Reference  10. 

To  obtain  the  necessary  input  tables  for  CMA,  the  user  selects  sets  of  values 
for  the  pressure  P,  the  dimensionless  gas  rate  B^,  and  the  dimensionless  char 
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rate  B^.  He  specifies  trie  elemental  composition  of  the  environment  gas,  the 
char,  and  the  pyrolysis  gas,  and  he  supplies  some  general  species  thermochemical 
data  for  all  molecules  he  wishes  considered  in  the  system.  Finally  he  specifies 
an  unequal  diffusion  coefficient  if  he  believes  unequal  diffusion  effects  will 
ba  important.  The  thermochemistry  code  then  computes  all  the  dependent  quanti¬ 
ties  of  interest  at  each  table  point  in  the  P  x  x  matrix  of  independent 

variable  values,  namely,  T  ,  Vz*  h.w,  and  h  ,  and  punches  this  information  out 

w  l**  iw  1  w  T 

on  punched  cards.  Similarly,  the  tables  of  Szie^iW  and  va^ues  f°r  the 
frozen  edge  gas  are  prepared  as  functions  of  P  and  T,  and  punched  out  on  cards. 
All  these  cards  form  part  of  the  card  input  deck  (see  Section  3.JL.10.3  below). 

2. 4. 3. 3  Simpler  Forms  of  the  Surface  Energy  Balance  Equation 

As  noted  in  Section  2. 4. 3.1  above,  for  equal  diffusion  coefficients 
the  Z*  driving  forces  reduce  to  the  simple  mass  fractions  K^.  If  in  addition 
to  equal  diffusion  the  user  specifies  that  Peuecjj  =  PeueCH'  t*ien  slnce 
£XiehjW  "  bew  and  £KiwhiW  =  hw  by  definition'  Equation  (9)  simplifies  to  the 
more  fiuailiar  form 

peueCH(Hr  ~  (1  +  B’)hw)  +  *chc  +  “ghg  +  awqrad  "  FloeTw  "  qcond  =  0  (12) 

T 

In  this  expression  hew  and  L2iehi.W  do  not  aPPear'  hence  the  frozen  edge  gas 
table  is  not  necessary  and  need  not  be  included  in  the  input  (see  Section 

3.1.10.3  below) . 

If  ablation  and  pyrolysis  do  not  occur.  Equation  (12)  assumes  a  still 
simpler  form 

peueCH (Hr  "  V  +  °wqrad  '  FlaeTw  "  qcond  =  0  (13) 

Section  3.1.10.3  below  will  make  it  clear  that  this  equation  can  be  handled  by 
the  CMA  code  through  the  mechanism  of  "non-ablating",,  surface  thermochemistry  •  ^ 
tables.  Equation  (13)  can  be  reduced  to  a  temperature  driving  force,  if  this 
is  desired 


h(Tr  -  V  +  awqrad  "  FloeTw  “  qcond  =  0 


(14) 


provided  that  the  user  enters  Tr  instead  of  Hr  and  h  instead  of  PeuecH  1°  the 
boundary  condition  time  tables  'xee  Section  3.1.9  below),  and  Tw  instead  of  hw 
in  the  surface  thermochemistry  tables  (see  Section  3.1.10.3  below). 
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2  4.3.4 


Automatic  Program  Adjustments  to  Input  Values  of  Time- 
Dependent  Convective  Heat  Transfer  Coefficient  PeueCH 

' 2.4. 3.4.1  General  remarks 

The  CMA  program  user  may  call  for  several  automatic  adjustments  to 
be  made  to  the  values  of  the  convective  transfer  coefficient  input  in  the 
functions-of-time  table.  This  section  describes  these  adjustments  and  the  order 
in  which  they  are  applied  in  the  program. 

One  adjustment,  pertinent  to  the  exposed  throat  ablative  materials  in 
rockets,  is  the  “radius  ratio  correction."  This  adjustment  accounts  for  the 
effect  of  pressure  decay  and  throat  radius  change  on  the  throat  transfer  coeffi¬ 
cient,  provided,  as  will  be  explained  below,  that  the  trends  of  the  effects  are 
given  by  the  Bartz  equation  (Reference  11) .  A  refinement  to  this  basic  correc¬ 
tion  includes  the  additional  effect  of  pressure  on  the  burning  rate  of  a  solid 
propellant.  A  second  refinement  on  the  throat  radius  ratio  correction  allows 
in  an  approximate  manner  for  the  effects  of  char  swelling. 

A  second  correction  allows  for  the  reduction  in  transfer  coefficients 
due  to  the  transpiration  or  blowing  effect  of  the  pyrolysis  gases  and  thermo- 
chemically  eroded  char  being  injected  into  the  boundary  layer. 

A  final  modification  distinguishes  the  convective  mass  transfer  coef¬ 
ficient  p  u  CM  from  the  heat  transfer  coefficient  p  u  C  by  referring  to  an  in- 

C  G  H  C  tl  » 

put  constant  ratio  Cjj/Cg- 

2. 4. 3. 4. 2  Radius  ratio  correction,  burning  rate  correction,  and 

char  swell  correction 

Radius  ratio  and  burning  rate  corrections 

If  the  user  is  analyzing  a  solid  propellant  rocket  throat  ablation 
problem,  and  calls  for  the  throat  radius  ratio  correction  according  to  the 
input  procedures  described  in  Section  3.1.10.3  below,  the  program,  will,  at 
every  step,  automatically  reduce  the  input  time  table  value  of  CH  to 

account  for  the  effects  of  the  increase  in  throat  radius  resulting  from  abla¬ 
tion.  The  correction  is  derived  by  noting  that  for  a  solid  rocket  motor 

pc  “  (^)  2/(1  -  n)  (15) 

where  Pc  -  chamber  pressure 

r*  -  nozzle  throat  radius 

n  -  burning  rate  exponent  of  the  solid  propellant  when  the 
burning  rate  r,  is  represented  by  r  =  bp£ 
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From  the  simplified  Bartz  equation 


«.« 


h_  * 


.*8.2 


(16) 


where  hc  is  the  convective  heat  transfer  coefficient, 
tions  gives 


h 


c 


cc 


1 


i.«-o.2n 

i-n 


Combining  these  two  equa- 


(17) 


For  n  =  0,  this  equation  reduces  to 


(18) 


and  would  represent  the  correction  in  the  convective  heat  transfer  coefficient 
for  a  constant  mass  flow  rate.  This  is  usually  termed  the  radius  ratio  correc¬ 
tion. 

The  CMA  program  refers  the  time  table  values  of  CH  to  the  initial 
condition.  This,  if  CHo  represents  the  input  heat  transfer  coefficient  at  a 
given  time,  the  program  will  reduce  the  input  value  according  to 


(19) 


if  the  radius  ratio  correction  has  been  called  for. 

Char  swell  correction 

It  has  been  determined  empirically  that  certain  phenolic  materials 
swell  or  thicken  during  or  following  charring.  The  amount  of  swell  has  in  some 
cases  (References  12  and  13)  appeared  to  be  roughly  proportional  to  the  thick-  . 
ness  of  the  charred  material: 


SswelT  KTchar 

where  K  is  a  proportionality  constant. 


(20) 
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The  char  thickness  (tc)  is  here  defined  as  the  char  depth  below  the 
original  surface  (6  )*  minus  the  sukrface  recession  predicted  in  the  absence  of 
swelling  (S  ) 


Sswell  V 


(21) 


Therefore,  if  swell  occurs,  the  correct  recession  would  be  given  by: 


S  =  S  -  S  .  =  S  -  K  (6  -  S) 

p  swell  p  c  p' 


(22) 


or 


S 


S  +  (S  -  6  )  K 
P  P  C 


(23) 


Furthermore,  surface  recession  rate  would  be  given  by: 

S  =  (1  +  K)S  -  K<5  (24) 

P  c 

For  a  given  value  of  K,  these  terms  are  computed  in  the  CMA  code  and  output  in 
an  auxiliary  output  block.  In  addition,  if  adjustments  in  heat  transfer  coeffi¬ 
cient  by  a  radius  ratio  correction  are  being  computed  as  is  often  the  case  with 
prediction  of  rocket  nozzle  throat  response,  the  corrected  (after  swell)  value 
of  recession  is  utilized  in  the  radius  ratio  correction  calculation.  Otherwise 
this  char  swell  formulation  is  not  used  in  uhe  thermal  response  calculations. 


2. 4. 3. 4. 3  Blowing  rate  correction 

The  CMA  program  will  automatically  reduce  the  transfer  coefficient 
to  account  for  the  familiar  blowing  effect  according  to  the  equation 


(25) 


where 


,  4  2>Mini 
"  peueCH, 


The  distance  6C  is  defined  as  the  distance  between  the  original  location  of  the 
heated  surface  and  the  present  location  of  a  line  of  density  pc  +  rc(pp  -  pc) , 
where  rc  is  a  user-selected  input  constant  (see  Section  3.1.6  below). 


and 


amount  of  material  injected  into  the  boundary  layer,  as 
discussed  below 

an  input  number,  discussed  below 


m.  . 
m3 


X 


peueCH, 


transfer  coefficient  as  input  and  reduced,  if  required, 
by  the  radius  ratio  and/or  burning  rate  effect  discussed 
in  the  section  immediately  above. 


Further  details  of  this  correction  scheme  will  be  discussed  below. 


Values  for  X 

With  X  =  0.5,  this  scheme  gives  the  "classical"  blowing  correction 
often  expressed  as  (References  14  and  15) 


tnd  +  mini/peueCH) 
-:/P  U  C. 


in]'  *e  e  H 


(27) 


which  is  useful  for  a  wide  range  of  problems. 

Other  values  of  X  allow  the  user  to  fit  blowing  correction  curves 
of  Ch/CHi  versus  or  to  account  for  special  effects,  in  the  few  cases 

where  these  are  known  with  confidence,  such  as  molecular  weight  effects*  or 
variable  property  effects.  In  view  of  the  uncertainties,  it  is  recommended 
that  X  -  0.5  be  used  for  laminar  flow.  A  value  X  =  0.4  appears  to  correlate 
constant  properties  turbulent  data  somewhat  better.  The  parameter  X  may  be 
made  a  function  of  time,  if  desired,  to  simulate  transition  effects.  A  value 
of  X  =  0  suppresses  the  blowing  correction. 


Options  for  n»inj 

The  symbol  m^nj  represents  the  rate  of  mass  injected  into  the  bound¬ 
ary  layer  which  is  effective  in  changing  velocity  profiles  and  reducing  diffu¬ 
sion  events.  The  CMA  program  allows  three  general  choices -for  m^: 

1.  All  mass  transfer  is  effective  in  reducing  blowing 


m.  . 
xnj 


=  mc  + 


m 


2.  If  part  of  mc  represents  failing  species,  and  this  failing  com¬ 
ponent  mfaj^  appears  on  the  surface  thermochemistry  cards  (an 

ACE  output  option) ,  then  these  fail  rates  can,  at  the  option 
0f  the  user  (see  Section  3.1.10.3  below)  be  read  in. 


These  are  known  to  be  important  in  some  simple  injection  problems  but  with 
many  chemical  reactions  taking  place  the  situation  is  somewhat  obscure. 
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stored,  and  deducted  as  appropriate  when  computing  injection 
rates  for  the  blowing  correction 


=  in  -  m,  +  rn 
erg 


3.  It  can  be  assumed  that  the  pyrolysis  gas  is  injected  into  the 
boundary  layer  through  cracks  and  fissures  rather  than  evenly 
over  the  heated  surface;  in  this  case  the  gas  injection  pre¬ 
sumably  plays  no  role  in  the  blowing  reduction,  in  which  case 


if  Option  2  above  is  being  used.  This  fissure  model  is  dis¬ 
cussed  in  more  detail  in  Section  2. 4. 3. 5  below. 


The  CMA  program  automatically  computes  the  injection  rate  nu 
according  to  the  first  option  above  unless  the  fail  option  has  been  invoked 
according  to  the  input  rules  discussed  in  Section  3.1.10.3  below  and  unless 
the  fissure  model  has  been  called  for  as  described\in  the  same  section. 


Blowing  correction  applied  to  C„ 

_ _1 

Note  that  the  blowing  correction  is  applied  to  the  value  of  CH  after 
it  has  been  adjusted  to  account  for  the  radius  ratio  (including  burning  rate 
and  char  swell)  effect. 


■  2. 4. 3. 4. 4  Computation  of  convective  mass  transfer  coefficient  p  u „C„ 

c  v  n 

After  PeuecH  has  been  computed  from  the  input  time  dependent  value 
PeueCH  and  corrected ,  if  called  for,  for  radius  ratio  and  mass  transfer  effects, 
a  mass  transfer  coefficient  peueCM  is  computed  by  multiplying  peueCM  by  an  in¬ 
put  constant  ratio  C^/Cy.  See  Section  3.1.10.3  below. 

2. 4. 3. 4. 5  Possible  future  additional  corrections 

The  automatic  adjustments  to  input  values  of  PeuecH0  presently  in¬ 
cluded  in  the  CMA  program  include  the  largest  and  most  important  corrections 
the  user  is  likely  to  desire.  A  small  temperature  correction  might  occasionally 
be  of  interest  but  has  not  yet  been  added  to  the  program.  A  molecular  weight 
effect  in  the  blowing  correction  must  for  the  present  be  accounted  for  by  ad¬ 
justing  X;  this  probably  will  yield  adequate  accuracy  until  more  is  known  about 
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(1)  the  desired  molecular  weight  effect,  and  (2)  the  molecular  weight  of  the 
pyrolysis  gas  component  of  the  injected  material. 


,  2. 4. 3. 5  Fissure  Model  Option  for  Surface  Energy  Balance 

Equation  (9),  which  is  repeated  here  for  reference, 
includes  the  energy  efforts  associated  with  injecting  the  pyrolysis 


peueCH(Hr  “ 


h„  ) 


peueCM 


E-t.  - 


Z* 

iw 


T 

,)h/  - 


B'h. 


m„h 
c  c 


m  h 

g  g 


wHrad 


-  FaewT;- 


^cond 


=  0 


(28) 


gas  rate  into  the  surface  from  below  with  enthalpy  hg  and  "blowing"  a  cor¬ 
responding  flux  away  from  the  surface  with  enthalpy  hw  (included  in  the  B’h^ 
term).  However,  some  ablative  materials  crack  or  fissure  during  heating  and 
sometimes  it  is  desired  to  explore  the  effects  of  assuming  that  the  pyrolysis 
gas  follows  these  fissures  to  the  surface  and  is  blown  out  through  the  boundary 
layer  without  participating  to  any  great  degree  with  the  boundary  layer  events. 

Specifically,  in  the  "fissure  model",  the  rah  injection  energy  term  in  the 

g  g 

surface  energy  balance  is  matched  by  an  equal  ra^h^  outflow  term  as  the  gas 

passes  across  the  surface  without  interacting.  By  the  same  token,  the  blowing 

term  B'h  now  becomes  (B1  -  B')h_,  since  the  B'  blowing  effect  has  been  isolated 
w  g  w  g 

and  identified  with  enthalpy  hg  rather  than  hw, 

* 

Thus,  the  "fissure  model"  surface  energy  balance  equation  is 


peueCH (Hr  * 


h  ) 


2UeCM  Z>ie  " 


T 

xwV  - 


*ch« 


in  h_ 
c  c 


+  aw^rad  FaewTw  ^cond  ® 

Suitable  programming  has  been  implemented  in  the  current  CMA  code  to  allow 
utilization  of  this  model  with  only  minimal  adjustment  of  input  data. 
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The  user  may  call  for  this  surface  energy  balance  type  by  punching  the  fissure 
model  flag  in  the  lead  card  for  the  surface  thermochemistry  tables  (see  Section 
3.1.10.3  below) .* 

When  using  this  model,  the  user  generally  makes  two  other  modifica¬ 
tions  to  the  usual  input.  First  to  harmonize  the  in-depth  calculations  with 
the  idea  of  the  pyrolysis  gas  rising  up  cracks  and  fissures,  pyrolysis  gas 
thermal  pick-up  in  the  char  is  suppressed  by  artifically  “flattening"  the  gas 
enthalpy-temperature  relationship  (see  Section  3.1.8  below)  for  temperatures 

above  the  pyrolysis  temperature  (i.e.,  Cn  =  0  in  this  temperature  range).  Sec- 

vg 

ond,  the  surface  thermochemistry  tables  are  prepared  for  l  =  0  so  as  to  exclude  a 
any  surface  thermochemical  effects  of  the  pyrolysis  gases.  Since  the  interpola¬ 
tion  routine  requires  two  values  for  a  given  pressure,  the  -  0  table  is 
duplicated  and  the  duplicate  is  assigned  a  second,  artifical  value. 

< 

2. 4. 3. 6  Accounting  For  Condensed  Phase  Removal  (Melting  or  "Failing") 

The  surface  energy  balance  Equation  (9)  does  not  show  an  energy  term 
corresponding  to  condensed  phase  removal  (as  in  thin  layer  melting  or  failing) 
but-  failing  is  detected  by  the  ACE  surface  state  program  and  the  associated 
energy  effect  is  automatically  accounted  for  in  the  surface  thermochemistry 
tables.  The  ACE  User's  Manual ,  Reference  13,  describes  the  melting  or  failing 
features  of  that  code  and  provides  detailed  instructions  for  implementing  them. 

The  EST3  code  cannot  account  for  failing. 

Or 


This  flag  also  causes  Bg  not  to  be  included  in  the  blowing  reduction  calcula¬ 
tion  for  the  transfer  coefficients  (see  Section  2.4.3. i  above). 
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SECTION  3 


USER’S  GUIujs 

This  part  of  this  report  provides 

•  Detailed  user  oriented  input  instructions  for  the  CMA  code 

•  An  explanation  of  the  program  output 

•  Miscellaneous  specific  information 

3.1  INPUT  FORMAT 

The  input  to  the  Charring-Material  Thermal  Response  and  Ablation 
Program,  Version  3,  can  conveniently  be  divided  into  ten  parts.  All  ten  por¬ 
tions  of  the  deck  must  be' present  for  most  runs.  These  individual  parts  will 
be  described  in  the  following  subsections. 

3.1.1  Title  and  Heading  Information 

The  first  three  cards  of  the  data  deck  are  used  to  transmit  title 
and  heading  information  to  the  output.  The  first  72  columns  of  each  of  these 
cards  may  be  used  for  the  title,  the  alphameric  information  in  columns  61 
through  72  of  the  third  card  being  used  as  a  page  heading  on  all  pages  after 
the  fitst. 

v 

3.1.2  Internal-Decomposition  Kinetic  Data 

These  three  cards  supply  -the  constants  used  in  Equation  (2).*  They 
are  supplied'0 according  to  the  following  format. 


Column 

Format 

Data 

Units 

1 

11 

Alphabetic  character  A,  B,  or  C  referring  to 
the  material  components.  These  must  be  on 
the  first,  second  and  third  cards  of  this  set 
in  order  f 

11-20 

F10.5 

Initial  density  of  component  i,  pQ^ 

lb/ft3 

21-30 

F10.5 

Residual  density  of  component  i,  pr 

i 

lb/ft2 

31-40 

F10. 3 

Pre-exponential  factor,  B^,  of  Equation  (2) 

sec  1 

41-50 

FlO .  5 

Density  factor  exponent, 

— 

* 

Reference  4  cites  relevant  data  for  many  materials  of  interest. 


Column 

Format 

Data 

Units 

51-60 

E10.3 

Activation  energy  factor  (Eai/°R) 

°R 

61-70 

F10.5 

Minimum  temperature  of  reaction  zone 

°R 

71-79 

— 

Not  read 

— 

80 

— 

(Third  card  only)  Number  (0  to  5)  of  decom- 

— 

posing  (pyrolyzing)  back-up  materials  in 
this  problem 


In  this  table,  the  minimum  temperature  of  the  reaction  zone  is  in¬ 
cluded  only  as  a  means  of  reducing  computatior al  time.  Thu3,  if  a  nodal  tem¬ 
perature  is  below  this  value,  the  program  will  bypass  the  fairly  complex  density 
calculation  procedure.  Also  for  the  case  of  a  non-decomposing  reinforcing  ma¬ 
terial,  the  value  of  this  temperature  can  be  set  to  a  ridiculously  high  value, 
thus  assuring  no  decomposition. 

3.1.3  Output  Interval  Specification  and  General  Program  Constants 

Two  cards  are  used  to  provide  the  program  with  the  values  of  certain 
general  constants,  to  establish  the  time  intervals  for  which  output  will  be  ob¬ 
tained,  and  to  specify  thermocouple  and  isotherm  output,  if  desired.  (If 
thermocouples  and/or  isotherms  are  called  for,  additional  cards  are  needed,  as 
described. ) 


First  Card 


Column 

Format 

Data 

Units 

1-2 

12 

Number  of  node lets  per  node  (blank  implies  10) 

— 

3-5 

13 

Total  number  of  nodes  (maximum  of  100) 

— 

6 

11 

A  non-zero  punch  in  this  column  will  cause 
thermocouple  and  isotherm  output,  normally 
listed  with  regular  output,  to  also  be 
punched  as  cards 

1 

CO 

12 

Number  of  thermocouples  for  which  output  ih 
desired.  (The  depths  of  these  thermocouples 
are  specified  on  additional  cards  of  this 
group.)  Punches  here  will  cause  thermocouple 
temperatures  as  a  function  of  time  to  be  out¬ 
put  at  the  end  of  the  regular  program  output. 

9-10 

12 

Number  of  isotherms  for  which  output  is  de¬ 
sired.  (The  temperatures  of  these  isotherms 
are  specified  on  additional  cards  of  this 
group.)  Punches  here  will  cause  isotherm 
depths  as  a  function  of  time  to  be  output  at 
the  end  of  the  regular  program  listing. 

_  3-2 


¥■ 


First  Card  (concluded) 


Column 

Format 

'  Data 

Units 

11-20 

F10.5 

Initial  value  o£  time 

sec 

21-30 

F10.5 

Final  value  of  time 

sec 

31-40 

F10.5 

Initial  output  time  interval 

sec 

41-50 

F10. 5 

/ 

Second  output  time  interval 

sec 

51-60 

F10.5 

Third  output  time  interval 

sec 

61-70 

F10.5 

Maximum  time  step  permitted  under  any  circum¬ 
stances.  If  unpfmched,  this  limit  will  be 
set  equal  to  5.0  seconds 

sec 

71-80 

F10.5 

Blowing-rate  parameter X 

f  _ 

Second 

Card 

Column 

Format 

Data 

Units 

1-10 

F10.5 

<* 

Time  of  transition  from  initial  to  second 
output  time  interval 

sec 

11-20 ; 

F10. 5 

Time  of  transition  from  second  to  third  out¬ 
put  time  interval 

sec 

21-30 

F10. 5 

i 

Minimum  thickness  of  last  ablator  node.  When 
because  of  surface  recession,  the  last  node 
of  the  ablation  material  shrinks  below  this 
value,  it  is  combined  with  the  adjacent  abla¬ 
tor  node.  A  value  equal  to  the  second  node's 
thickness  is  usually  appropfiate 

in 

31-40 

41-50 

51-60 

F10.5  - 
*  F10.5 

F10. 5 

Heat  of  formation  of  the  virgin  plastic,  the 
char,  and  the  pyrolysis  gas,  respectively. 

These  are  evaluated  at  the  datum  temperature 
given  in  columns  71-80 

Btu/lb 

61-70 

F10.5 

Volume  fraction  or  mass  fraction  of  the  vir¬ 
gin  plastic  which  is  occupied  by  resin.  If  ■  * 
mass  fraction,  it  is  input  as  a  negative 
number  ^  . 

71*80 

F10.5 

Datum  temperature  for  heats  of  formation  given 
in  columns  31-50.  For  Option  1  computations, 
this  must  be  the  same  temperature  as  used  in 
the  EST  program  for  generating  the  surface 
thermochemistry  data  deck. 

°R 

3-3 


Additional  Cards  as  Required 


Column  Format 

1-80  8 (F10. 5) 


If  thermocouple  and/or  isotherm  output  have 
been  called  for,  these  cards  specify  thermo¬ 
couple  depths  in  inches  below  the  original 
surface  and  isotherm  temperatures  in  degrees 
Rankine.  Isotherm  temperatures  must  begin  on 
a  new  card .  The  total  number  of  thermocouples 
and  isotherms  may  not  exceed  20.  If  no  ther¬ 
mocouples  or  isotherms  are  called  for,  these 
card(s)  are  omitted. 


Units 

in 

and/or 


3.1.4 


Back-Up  Material  Decomposition  Data  (Optional) 


If  the  problem  includes  decomposing  (pyrolyzing)  back-up  materials, 
a  four  card  set  of  pyrolysis  data  must  be  included  for  each  such  decomposing 
back-up  material.  Each  four  card  set  has  two  subgroups: 

o  First  Subgroup: 

Back-Up  Material  Decomposition  Kinetic  Data  (3  Cards) 


Three  cards  in  the  format  of  the  three  kinetic  data  cards  for  the 
main  material,  except  that  column  80  of  Card  3  is  not  read.  Groups 
4a  (three  cards)  and  4b  (one  card)  are  placed  together  as  a  four- 
card  unit.  There  will  be  one  such  four-card  unit  for  each  decom¬ 
posing  back-up  material. 

Second  Subgroup: 

Enthalpies  of  Formation,  Resin  Volume,  Or  mass  Fraction,  En 
Enthalpy  Reference  Temperature  (1  Card) 

This  card  follows  the  same  format  as  the  second  card  of  Section^ 
3.1.3  except1- that  general  data  on  that  card  not  pertinent  to  the 
back-up  material  is  not  needed  and  is  not  read; 


Column 


31-40 

41-50 


51-60 


61-70 


Format 

ad  lib 

F10 . 5 
F10.5 


ad  lib 


F10 . 5 


Not  reaa 

Heats  of  formation  of  virgin  plastic  and  char, 
respectively,  evaluated  at  the  datum  tempera¬ 
ture  given  in  columns  71-80 

Not  read 

\ 

Volume  fraction  or  mass  fraction  of  the  virgin 
plastic  back-up  which  is  occupied  by  resin.  If 
mass  fraction, .-i,t„j.s  input  as  a  negative 
number 


Units 


Btu/lb 


71-80 


F10 . 5 


Datum  temperature  for  heats  of  formation. 


3-4 


3.1.5 


Nodal  Data 

A  set  of  cards  equal  in  number  to  the  number  of  nodes  is  used  to 
provide  certain  information  with  regard  to  the  initial  state  of  the  nodes, 
their  thickness,  and  the  contact  resistances  between  them.  In  addition,  the 
first  card  of  this  set  is  used  to  specify  the  geometric  nature  of  the  exposed 
surface,  that  is  whether  it  is  a  flat  plate  or  an  internal  or  external  radius. 
The  format  for  these  cards  is  shown  in  the  following  table. 


Column 

Format 

Data 

Units 

1-2 

12 

Material  number,  1  for  main  material  virgin 
plastic,  2  for  main  material  char,  0  for  main 
material  pyrolysis  zone  (can  not  be  input) , 

22-31  for  charring  back-up  materials  (see  dis¬ 
cussion  below) ,  and  3-10  for  non-charring  back¬ 
up  materials 

3-12 

E10.0 

Initial  temperature  of  the  node 

°R 

13-22  £10 .0  Initial  cross-sectional  area  of  the  node  in 

tiny  units  or  normalized  on  any  convenient 
reference  such  as  the  surface  area.  This 
data  is  required  only  for  "general"  geome¬ 
tries;  for  planar,  cylindrical,  and  spherical 
geometries,  and  for  shapes  for  which  the  area 
varies  as  a  power  of  the  radius  the  specifica¬ 
tion  of  individual  nodal  areas  is  not  required 
and  these  columns  should  be  left  blank 

23-32  E10.0  Initial  thickness  of  the  node  (For  suggestions 

v  on  the  selection  of  nodal  thicknesses,  see 
Appendix  A.) 

33-42  E10.0  First  card  only*  0  for  planar  surface;  posi- 

tive  value  of  Internal  radius  or  negative  of 
value  of  external  radius  for  all  non-planar 
axisymmetric  geometries  (including  spheres 
and  shapes  with  A  ~  ra)  in 

Second  card  only.  Exponent  on  radius  giving 
nodal  cross-sectional  area  variation  for 
spherical  and  other  "regular"  geometries. 

Sphere  requires  2.0.  No  entry  needed  for 

planar  and  cylindrical  geometries.  Must  be 

blank  is  "general"  geometry  option  is  used  in 

43-52  E10.0  Contact  resistance  between  this  node  and  the  f tzsec-°R/Btu 

next  node. 

53-80  38X  NOT  USED 


Arbitrary 
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The  number  of  nodes  may  not.  exceed  100.  Note  that  for  a  restart  cal¬ 
culation,  the  nodes  that  make  up  the  precharred  depth  in  the  main  material  are 
called  out  as  material  2. 

In  the  list  of  nodes,  any  decomposing  back-up  material  nodes  will 
appear  between  the  last  node  of  the  main  ablating  material  and  the  first  node 
of  non-charring  back-up  (if  any) .  Material  numbers  of  the  charring  back-ups 
must  be  assigned  as  follows: 


Virgin 

Char 

1st  charring  back-up 

22 

23 

2nd  charring  back-up 

24 

25 

3rd  charring  back-up 

26 

27 

4th  charring  back-up 

28 

29 

5th  charring  back-up 

30 

31 

The  rules  here  differ  from  those  for  non-charring  back-ups  in  that  one  charring 
back-up  must  be  distinguished  from  another  if  it  has  a  different  location.  For 
example,  the  fourth  back-up  may  happen  to  be  identical  to  the  first,  but  it 
still  receives  a  different  material  number. 

3.1.6  Back  Wall  Heat  Transfer  Conditions;  Char  and  Pyrolysis 

Zone  Criteria 

Heat  transfer  at  the  back  wall  is  characterized  by  a  convective  heat 
transfer  coefficient,  an  emissivity,  and  the  temperature  of  a  "reservoir”  to 
which  heat  transfer  takes  place  from  the  back  wall.  This  card  allows  the  spec¬ 
ification  of  these  three  quantities.  Blanks  serve  to  specify  an  insulated 
back  wall. 

Two  additional  entries  define  char  edge  density  and  pyrolysis  zone 
edge  density  according  to  the  following  definitions: 

pchar  =  pc  +  rc(pP  -  pc> 
edge 

p  i  .  =  p  +  r  (p  -  p  ) 

^pyrolysis  c  p  p  c 

edge 

The  user  inputs  values  for  rc  and  r^.  Typical  values  are  rc  =  0.02  (defining 

the  char  edge  as  occurring  where  the  density  is  pc  +  0.02 (Pp  -  pc> )  and  r^  = 

0.98.  If  no  entries  are  made  for  r  and  r  .  0.02  and  0.98  will  be  assumed. 

c  p 
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For  char  swell  corrections,  a  value  of  rc  =  0.5  is  suggested  (see  Section 
2.4. 3.4  above) . 


Column 

Format 

Data 

Units 

1-10 

F10.5 

Back  wall 

convective  coefficient 

Btu/sec-ft2  °R 

11-20 

F10.5 

Back  wall 

emissivity 

— 

21-30 

F10.5 

Reservoir 

temperature 

°R 

31-40 

F10.5 

Char  zone 

criterion  r 

c 

— 

41-50 

F10.5 

Pyrolysis 

zone  criterion  r 

P 

— 

60 

11 

One  punch  calls  for  output  of  thermal  con¬ 
ductivity  in  place  of  enthalpy  in  standard 
output  block 

3.1.7 

Material 

Property  Tables  and  f -Function  Tables 

In  these  tables  are  presented  the  requisite  thermodynamic  data  for 
the  main  material  virgin  plastic  (material  no.  1) ,  the  main  material  char  (ma¬ 
terial  no.  2),  any  charring  back-up  materials,  and  the  various  non-charring 
back-up  materials  (materials  3  through  10,  as  required)  specified  in  the  nodal 
data  table.  These  data  are  input  as  functions  of  temperature  according  to  the 
following  format. 


Column 

Format 

Data 

1-2 

12 

Flag,  nominally  zero,  +1  marks  terminal  card 
of  last  material  property  table,  -1  marks 
terminal  card  of  other  intermediate  material 
property  tables. 

3-12 

F10.5 

Temperature  (independent  variable) 

13-22 

F10 . 5 

Specific  heat 

23-32 

F10. 8 

Thermal  conductivity 

33-42 

F10 .5 

Emissivity 

Units 


o 


R 


Btu/lb-°R 

Btu/ft-sec-°R 


Separate  tables  must  appear  for  the  main  material  virgin  plastic  and 
the  main  material  char,  in  that  order.  If  the  main  material  does  not  char,  a 
"char"  table  must  still  be  included.  If  charring  back-up  materials  have  been 
included  in  the  problem,  then  paired  tables  must  be  given  for  the  virgin  plastic 
and  char  properties  of  all  such  back-ups,  with  the  table  pairs  appearing  in  the 
same  order  with  which  such  back-ups  appear  in  the  nodal  list  (see  Section  3.1.5 
above.  If  a  charring  back-up  is  duplicated,  duplicate  tables  must  be  included 
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in  the  proper  order.  Tables  for  decomposing  back-ups  follow  the  two  tables  for 
the  main  charring  material  and  precede  any  tables  for  non-charring  back-ups. 

Any  necessary  property  tables  for  non-charring  back-up  materials 
follow  these  charring  material  tables.  Preceding  each  such  table,  a  single 
lead  card"  is  used  to  specify  the  back-up  material  identification  number  and 
density.  These  data  are  given  in  the  first  ten  Columns  with  the  following 
format. 


Column 

Format 

Data 

Units 

1-2 

12 

Material  identification  number  (ranging  from 
3-10,  inclusive) 

— 

3-10 

F8.4 

Material  density 

lb/ft3 

If  the  emissivity  for  materials  1  and  2  is  input  as  zero,  it  is  auto¬ 
matically  presumed  to  be  1.0  for  any  Option  3  computations.  (The  user  may  de¬ 
sire  to  input  a  zero  emissivity  for  radiation  equilibrium  problems  under  Option 
1.) 

The  number  of  temperature  points  in  each  table  may  not  exceed  30  or 
may  not  be  less  than  2.  The  tables  must  be  ordered  on  either  ascending  or  de¬ 
scending  temperatures. 

For  all  charring  materials,  material  properties  for  partially  de¬ 
graded  or  charred  material  are  formed  from  weighted  averages  of  the  plastic 
and  char  properties  for  the  relevant  nodal  temperatures  according  to  the  mix¬ 
ture  rules  presented  in  Section  2.1.2  above.  If  the  user  wishes  to  employ  the 
so  called  "f-functions"  for  this  purpose  (see  Equation  (11)  of  Section  2.1.2 
above,  the  f ^  (x)  and  f 2 (x)  values  are  input  along  with  the  other  material  prop¬ 
erties  tables  described  in  this  section.  Each  f-function  table  follows  the 
following  special  rules: 

•  An  f-function  table  is  flagged  by  a  lead  card  of  the  following 
format 


Column 


1-2 


3-10 


Format  Data 

12  Identification  flag  which  must  be  between 

3  and  10  and  j^hich  must  not  duplicate  the 
identification Npumber  of  a  previously  read 
non-charring  back-up  property  table 

F8.4  Must  be  blank  or  zero 


Units 
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Column 


Format 


Data 


Units 


11-16  611  Material  assignment  numbers  for  the  follow¬ 

ing  f-f unction  table.  If  this  table  is 
assigned  to  the  surface  material,  column  11 
must  be  blank.  If  the  same  f-function  table 
is  to  be  assigned  to  one  or  more  charring 
back-up  materials,  these  assignments  follow 
in  columns  12-16  in  any  order  with  .  imbedded 
or  following  blanks  being  ignored.  Note 
that  for  the  material  assignments  in  columns 
11-16  of  the  lead  card,  the  surface  material 
is  designated  as  zero  (or  blank)  and  the 
charring  back-up  materials  are  designated  1 
through  5  in  the  same  order  as  they  appear 
in  the  nodal  data  cards.  If  column  11  is 
not  blank,  the  f-function  table  will  not  be 
used  for  the  surface  material. 


The  actual  f-function  table  follows  the  lead  card  and 


has 

i  the  following  format 

Column 

Format 

Data 

Units 

1-2 

12 

\ 

Flag,  as 'above,  for  last  card  of  table 

— 

3-12 

F10. 5 

x  (independent  variable) 

— 

13-22 

F10. 5 

fx(x) 

23-32 

F10. 5 

f2(x) 

— 

* 

•  The 

f-function  tables  may  be  distributed  anywhere 

among  the 

non- 

charring  back-up  properties  tables. 

•  If  any  charring  material,  main  or  back-up,  is  not 

assigned  an 

f-function  table,  the  thermal  conductivity  of  partially  charred 
substance  for  this  charring  material  will  follow  the  simple 
linear  form  of  Equation  (14). 


3.1.8  Pyrolysis  Gas  Enthalpy  Table 

This  table  specifies  the  variation  of  pyrolysis  gas  enthalpy  with 
temperature .  These  values  are  added  to  the  heat  of  formation  of  the  pyrolysis 
gas  as  specified  on  the  program  constant  cards  (see  Section  3.1.3  above).  This 
table  is  presented  on  card  pairs,  the  first  of  which  presents  a  set  of  tempera¬ 
ture  (up  to  8),  the  second  containing  the  corresponding  enthalpy  values.  The 
first  column  of  the  temperature  card  is  used  as  a  flag,  nominally  blank  but  con¬ 
taining  an  integer  to  mark  the  last  pair  of  cards.  The  format  for  this  card 
pair  is  given  as  (II,  F9.5,  7F10 . 5/8F10 . 5) ,  implying  a  basic  field  length  of 
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ten  columns  for  both  temperature  and  enthalpy,  with  the  exception  of  the  first 
temperature  which  is  restricted  to  columns  2  through  10.  The  temperatures 
should  be  given  in  ascending  sequence  with  the  exception  of  the  last  card  pair 
which  may  have  from  one  to  eight  entries).  The  total  number  of  temperature 
points  in  the  table  may  not  exceed  30  and  may  not  be  less  than  2. 

Note  that  no  pressure  dependence  of  pyrolysis  gas  enthalpy  may  be 
accounted  for.  Also  note  that  only  one  pyrolysis  gas  is  allowed,  and  hence  the 
gas  given  off  by  any  charring  back-ups  must,  be  the  same  as  the  pyrolysis  gas  of 
the  main  material,  as  would  be  the  case,  for  example,  for  carbon  phenolic  backed 
by  silica  phenolic. 


■* 


A 


3.1.9  Surface  Time-Dependent  Boundary  Conditions 

The  table  of  time-dependent  boundary  conditions  is  used  to  specify 
alternative  sets  of  time-varying  dependent  variables  for  the  various  options. 
The  table  must  be  ordered  with  £he  time  values  increasing.  -  ' 

The  table  may  be  thought  of  as  consisting  of  a  number  of  sub-tables 
each*  sub- table  representing  one  option.  The  switch  from  dne  option  to  another 
requires  a  repeated  time  entry,  the  first  card  representing  the  last  entry  of 
the  earlier  table  and  the  second  card  representing  the  start  c£  the  next  table. 

The  total  number  of  time  values  in  the  table  is  limited  to  30.  Each 
sub-table  must  have  at  least  2  entries,  hence  the  number  of  sub-tables  cannot 
exceed  16.  (The  most  common  problem  has  only  2,  representing  an  Option  1  or 
Option  2  calculation  followed  by  cooldown.  Option  3.) 

The  format  for  the  time-tables  is  as  follows: 


Column 

Format 

Data 

1 

11 

Flag,  nominally  blank,  punched  to  indicate 
the  last  card  of  the  time  table.  In  any 
Option  1  calculations  are  to  be  done,  this 
flag  must  be  a  one  (which  calls  for  reading 
the  surface  equilibrium  tables) .  If  no  Op¬ 
tion  1  calculations  are  to  be  done,  this 
flag  must  be  a  two. 

2-10 

F9 . 5 

Time  (independent  variable) 

11-20. 

F10.5 

Option  1:  /Recovery  enthalpy,  relative  to 
the  same  chemical  base  state  as  used  with 
the  heats  of  formation. 

Option  2:  Surface  temperature 

Option  3:  Radiation  view  factor  i 

Units 


Lee 

Btu/lb 

°R 


N. 

3-10 


* 


I 

I 


Columns 

Format 

/- 

Data 

Units 

21-30 

F10 .5 

Option 

Is 

Radiant  energy  flux  to  the  surface  * 

Btu/ft2-sec 

Option 

2: 

Surface-recession  rate 

mils/sec 

Option 

3: 

Radiant  energy  flux  to  the  surface 

Btu/ft2-sec 

31-40 

F10 . 5 

Option 

1: 

Heat-transfer  coefficient 

lb/ft2-sec 

Option 

2: 

Blank 

Option 

3: 

Must  be  blank 

41-50 

F10.5 

Option 

1: 

Pressure 

atm 

51-60 

F10 .5 

Option 

1: 

Blowing  reduction  parameter  (if  a, 
function  of  time;  blank  entries  ; 
will  be  filled  by  constant  value 
entered  in  column  71-80  of  the 
first  card  of  Section  3.1.3  above). 

3.1.10  Surface  Thermochemistry  Daca  (Option  1  Only) 

3 .1 .10.1  Introduction 

As  discussed  in  Section  2.4.3  above,  problems  involving  use  of  the 
surface  thermochemistry  option  (Option  1)  require  the  input  of  an  arry  of  sur¬ 
face  ^mass  and  energy  data  particular  to  the  option  and  the  material  Deing  ana¬ 
lyzed.  These  data  include  the  specification  of  the  ratio  of  mass  transfer  co¬ 
efficient  to  heat  transfer  coefficient,  the  radiation  view  factor,  and  surface 
thermochemical  data.  The  cards  containing  this  input  are  described  below. 

Problems  not  involving  Option  1  calculations  do  not  need  the  surface  thermo- 
1  /  • 
chemistry  data  deck. 

3.1.10.2  Ratio  of  Mass  to  Beat  Transfer  Coefficient,  View  Factor, 

'  Surface  Option  Flags 

A  single  card  serves  to  specify  the  ratio  cf  the  convective  mass 
transfer1 coefficient  to  the  convective  heat  transfer  coefficient  (CM/CH) •  the 
surface  radiation . view  factor,  an,d  a  number  of  control  flags.  (The  view  factor 
on  this  card  will  only  be  used  for  Option  1  computations.  For  Option  3  (cool¬ 
down)  computations,  a  separate  Option  3  view  factor  is  input  as  a  function  of 
time  in  the  functions-of-tkme  tables.  In  Option  I  problems  with  "radiation 
equilibrium* ,  the  radiant  energy  flux  to  the  surface  may  be  input  as  zero,  re¬ 
quiring  then  that  the  Option  1  view  factor  be  input  as  zero.) 


/• 

<■ 


Column 

Format 

Data 

Units 

’  1-10 

FlO.O 

CM/CH 

- A 

11-20 

FlO.O 

Option  1  view  factor 

— 

21-29 

F9.0 

Burning  rate  exponent  n  used  in  throat  prob- 

— 

lems  to  adjust  input  convective  transfer  co¬ 
efficient  PeueCH  according  to 

<i 


30 


31-39 

40 


50 


/ 


51-59 

60 


CH^  is  then  adjusted  to  account  for  blowing 

effects.  Applies  to  inner  radius  heating 
problems  (nozzles)  only.  See  Section 
/'2.4.3.4  above.  . 

/»  /  t 

II  One  punch  calls  for  radius  ratio  correction 

on  ;  nput  convection  transfer  coefficient 
peuecH0  according  to  Equation  (30)  above... 

Note:  Input  of  n  >  0  colls  for  Equation 

(30)  and  will  -  override  i.'r..k  in  this 
column.  'Applies  to  inner  radios  heating 
,  .  problems  (nozzles)  only. 


9X 


‘  Blank 


II 


II 


Pne-punch  calls  for  reuse  of  previously 
input  surface  tables.  No  more  input  is  r 
read.  Note:  Ratio  of  Cm/Cjj  (columns  1-10) 
must  be  the  same  for  thi3  calculation  as 
for  previous  one  if  this  option  is  to  be 
used.  P 

x. 

One-punch  reads  Bj  in  "new  format"  surface 
thermochemistry  tables  (as  discussed  below) 
and  uses  B'  -  B£  in  blowing  correction  to 
PeueCH  or  B'  -  Bg  —  B£  if  fissure  mpdOl  is 

being  used,  se~.  column  60  below)  ,  no  punch 
read’s  Told  format"  surface  tables  or  ignores 

Bf 


in  "new  format"  tables  and  uses  B*  in 


.blowing  correction  to 
fissure  model) .  ‘ 


peueCH 


(or  B'  '-  B^  xn 


9X 

II 


Blank 

One-punch  invokes  "fissure  nuj>del":  exludes 
Bg  from  blowing  correction  on’  peUeC^  and 

f>-om  surface  energy  balance;  see  Section 
2 . 4 . 3 . 5  above .  * 

Char  swell’  proportionally  constant  K  in 
Equation  (21),  section/2.4.3.4.2  '  -  - 

*  .  ’V 


/ 


3-12 


>  * 


61-70 


FlO.O 


Il 


3.1.10.3  Surface  Thermochemistry  Table 
3.1.10.3.1  Introduction 

This  table-  supplies  the  necessary  input  data  for  the  surface  energy 
balance  computations  in  Option  1.  (This  energy  balance  is  discussed  in  Section 

2.4.3  above.)  , 


Most  commonly  the  deck  of  cards  which  make  up  the  surface  thermochem¬ 
istry  table  is  generated  by  the  Equilibrium  Surface.  Thermochemistry  (EST)  pro¬ 
gram  or  the  Aerotherm  Chemical  Equilibrium  (ACE)  program.  The  user's  manual  for 
these  programs  (References  8  and  9  )  describe  this  table  in  complete  detail. 

On  occasion,  the  user  may  desire  to  construct  his  own  surface  thermochemistry 
table,  and  so  the  following  sections  include  brief  descriptions  of  the  organiza¬ 
tion  and  format  of  these  tables.  ;  The  main  emphasis,  however,  will  be  on  the 

•  •  .  *  !  ,  ? 
communication  between  *  therr.tochenis  try  programs  and  the  Charring  Material 

Ablation  program,  sin'ce  this  is  of  the  most  general  interest. 

■>  *  C3" 

3.1.10.3.2  Edge:  enthalpy  data  . 


Equation  (9)  of  Section  2.4.3  shows  that  if  diffusion  coefficients 
are  not  equal  or  if  the 'ratio  C^/C^  is  not  unity,  then -the  surface  ene’ISgy  bal¬ 
ance  requires  da*  about  the  edge  gases  of.  -the  boundary  layer.  These  dita  ape 


provided ^iti ^a  spe -  al  "edge  table".  The  independent  variables  for  this  table 
,are  pressure  ano  temperature . ,  Dependent  variables  are  he  \nd  the  sum  IZtgh^w. 


The  - .  .  enthalpy  data  are  entered  on  the  cards  as  follows: 


c 

Column 

Format 

) 

/ 

/ 

Data 

Units-  4 

1-8- 

F8.5 

Pressure. 

atm 

r 

9-24 

16X 

Blank 

25-33  ' 

F9.4 

Temperature 

°K 

(°R  if  Negative, 
in  which  case 
enthalpies  be¬ 
low  are  Btu/lb) 

.34-38 

39-47  ' 

*»* 

F5.3 

F9.3 

Unequal  diffusion  exponent 

T 

Summation  Tzt  h,w 
t  le  l 

cal/gr 

(Btu/lb  if  tem¬ 
perature  is  en- 
entered  with 
minus  sign) 

48-56 

F9.3 

<> 

\ 

cal/gr 

(Btu/lb  if  tem- 
•  perature  is 
entered  with 
minus  sign) 

57-58  12 


-1  (flag  signifying  that  this  card  is 
part  of  the  edge  gas  table) 


t 


3-13 


Column 

Format 

Data 

Units 

59-66 

2A4 

Unused,;  but  may  contain 
mation 

alphameric  infor- 

- 

67-78 

2A6 

Problem  identification 

(not  read) 

— 

'79-80 

2X 

UNUSED 

Note  that  although  the  EST  and  ACE  programs  will  provide  a  data  deck  usinij  °K 
and  cal/gr ,  in  those  rare  cases  in  which  a  user  wishes  to  supply  his  own  deck  / 
and  prefers  to  work  in  °R  and  Btu/lb,  he  may  do  so  simply  by  introducing  a  ./ 
minus  sign  as  a  flag  in  front  of  the  temperature  entries. 

The  table  length  is  limited  to  5  pressure  sets  (it  may  have  only  1 
pressure  set)  with  noc  more  them  30  nor  less  them  3  temperature  entries  in  each 
sjet.  The  series  of  temperature  values  may  be  different  for  each  pressure  set. 
The  table  has  been  organized  as1' a  series  of  sections,  each  representing  one 
pressure  and  each  preceding  the  corresponding  pressure  group  of  the  surface 
thermochemistry  deck  as  described  below.  The  temperature  entries  within  each 
section  must  be  ordered,  either  ascending  or  descending.  Similarly,  the  pres- 
•sures  must  be  ordered  either  ascending  or  descending.  (Deck  generated  by  the 
EST  and  ACE  programs  will  have  been  automatically  ordered  properly.)  The 
sketcn  shown  on  the  following  page  demonstrates  t!he  thermochemical  data  tables 
make-up .  f 

3.1.10-3.3  Surface  thermochemistry  table 
“  This  table  is  comprised  of  a  se:ries  of  sections  of  up  to  29  cards 

each.  The  sections  represent  one  pressure  itnd  one  pyrolyses  gas  rate.  Each 
section  consists  of  two  subsections.  The  first  contains  the  ablating  cases; 
here  the  char  ablation  rate  is  the  third  independent  variable  and  the  surface 

I  .C 

temperature  is  a  dependent  variable.  The  second  subsection  represents  surface 
temperatures  too  low  for  ablation;  in  this  subsection  surface  temperature  is  the 
independent  variable.  ° 

Thus  one  table  has  three  independent  variables:  pressure,  pyrolysis 
gas  rate,  and  either  surface  temperature  or  ablation  rate,  depending  on  whether 
the  surface  temperature  is  high  enough  for „ ablation. 

The  table  has  either  two  or  three  dependent  variables,  according  to 

whether  the  surface  temperature  is  high  enough  for  ablation.  Two  dependent 
*  T 

variables  always  present  are  the  summation  I2isvhjW  an^  hw'  t*ie  enthalpy  of  the 

wall  ijases.  The  third  dependent  variable  is  the  surface  temperature,  but  it  is 

dependent  only  in  those  cases  for  which  the  surface  temperature  is  high  enough 

for  ablation.  (Otherwise,  the  surface  temperature  functions  as  an  independent 

variable . ) 


The  EST  and  ACE  programs  generate  separate  groups  for  each  pressure, 
one  at  a  time.  These  groups  must  be  ordered  on  pressure  (either  ascending  or 
descending)  by  the  user  to  make  up  the  surface  thermochemistry  deck.  With  each 
pressure  group  the  pyrolysis  gas. rates  Py  will  be  ordered  in  descending  order. 
Within  each  gas  rate  section;’’  non-zero  char  rate  entries  will  be  grouped  ahead 
of  the  zerb  char  rate  entries.  The  non-zero  char  rate  entries  will  not  be  or¬ 
dered  in  any  particular  way  on  the  char  rates;  any  necessary  ordering  is  made 
automatically  by  the  Charring  Material  Ablation  program  as  it  reads  the  data. 
The  zero  char  rate  entries  are  ordered  with  descending  temperatures.* 


(Users  providing  their  own  thermochemistry  decks  must  ensure  that  the 
pressures  and  gas  rates  are  ordered,  but  the  ordering  may  be  either  ascending 
or  descending  in  each  case.  Within  each  gas  rate  section,  non-zero  char  rate 
entries  must  precede  the  zero  char  rate  entries  but  need  not  be  ordered.  Char 
rates  may  not  be  duplicated  in  a  given  table.  Zero  char  rate  entries  follow 
and  must  be  ordered  with  descending  temperatures.  These  cards  are  identified 
as  zero  char  rate  cards  by  a  flag  in  columns  57-58,  as  described  in  the  format 
specification  below.) 


The  number  of  pressure  groups  may  not  exceed  5  (and  may  be  only  1) ; 
the  number  of  pyrolysis  gas  rates  in  each  pressure  group  may  not  exceed  20  and 
may  not  be  less  than  2.  The  sequence  of  gas  rate  values  must  be  the  same  in 
the  different  pressure  sections/  Within  each  gas  rate  group  the  number  of  char 
rate  entries,  including  the  zero  char  rate  (independent  surface  temperature) 
cards,  may  not  exceed  30  and  may  not  be  less  than  2.  The  series  of  char  values 
may  be  unique  for  each  section.  j 

(The  °R-Btu/lb  option  described  for  the  edge  tables  in  Section 
3.1.10.3.2  may  be  used  for  these  tables  also.) 


The  card  format  for  the  surface  thermochemistry  data  is  as  follows: 

i 


Column 

Format 

Data 

1-8 

F8.5 

Pressure 

9-16 

F8.5 

Gas  rate  mg/PeueCM  = 

17-24 

F8.5 

Char  rate  m  /p  u  Cu  =  B' 
c  e  e  M  c 

25-33 

F9.4 

Surface  temperature 

Unit 

atm 


°K 

{  R  if  negative 
in  which  case 
enthalpies  be¬ 
low  are  Btu/lb) 


Limitations  in  the  EST  and  ACE  chemistry  routines  sometimes  require  that  thesi 
"zero  char  rates"  not  be  zero,  but  some  small  number,  for  example,  0.0001. 
This  causes  no  difficulty  in  CI>A. 


3-16 


/ 


Column ^  Format 
34-38  F5.3 

39-47  F9.3 


Data 

Unequal  diffusion  exponent 
T 

Summation  IZ*W^W 


48-56  F9.3 


Enthalpy  of  wall  gases  hw 


57-58  12 


59-60  2X 


0  for  assigned-temperature  entries  in 
the  thermochemistry  program  (no  abla¬ 
tion)  ;  >  0  for  standard  surface  thermo¬ 
chemistry  with  ablation  (temperature  is 
dependent) 

Blank 


61-66  A6 


67-78  2A6 

79-30  12 


Chemical  symbol  of  surface  species 
(EST  program  prints  such  symbols 
arranged  alphabetically  and  truncated 
from  right  end  if  necessary.) 

Problem  identification  (not  read) 
UNUSED 


Units 


cal/gr 

(Btu/lb  if  tem¬ 
perature  is  - 
entered  with 
minus  sign) 

cal/gr 

(Btu/lb  if  tem¬ 
perature  is 
entered  with 
minus  sign) 


The  ACE  program  can  account,  as  one  of  its  options,  for  mechanical 
removal  or  failing  of  candidate  surface  species  at  defined  temperatures .  The 
energy  associated  with  this  failing  cam  be  (and  is)  accounted  for  by  the  ACE- 
CM&  program  combination;  this  is  accomplished  by  ACE  by  suitably  modifying  the 

i 

hw  output  term  on  the  punched  card  thermochemistry  table  (thus.hw  does  not 
represent  the  enthalpy  of  the  wall  gases  if  failing  is  occurring) . 

The  actual  value  of  the  dimensionless  failing  rate  B£  =  ®f/peueCM 
may  or  may  not  appear  on  the  punched  surface  thermochemistry  card  tables,  de¬ 
pending  on  whether  or  not  the  ACE  user  has  called  for  the  appropriate  card  out¬ 
put  format.  This  "new  format"  differs  from  that  cited  above  in  columns  61-80: 

( 


Alternate  Card  Format  -  Surface  Tables 


Column 

Format 

Data 

Units 

1-60 

As  described  above 

61-68  , 

2A4 

Chemical  symbol  of  surface  species 

— 

69-78 

E10.3 

Bf  *  VpeueCM 

— 

79-80 

2X 

Blank 

— 

3-17 


The  CMA  user  may  call  for  the  CMA  program  to, read  these  B£  entries 
and  to  allow  for  them  in  the  blowing  correction  to  Chq  by  punching  a  one  in 
column  50  in  the  lead  card  for  the  surface  thermochemistry  tables,  as  described 
in  Section  3.1.10.2  above. 

3.1.10.3.4  Termination  card 

The  surface  thermochemistry  data  deck  must  be  terminated  by  a 
single  blank  card.  Output  decks  of  the  EST  and  ACE  programs  may  not  have  such 
a  card,  in  which  case  the  user  must  supply  *it. 


3.2.1 


PROGRAM  OUTPUT 

Input  Data  (Except  Surface  Thermochemistry  Tables) 


Program  output  begins  with  an  output  of  the  input  title  and  heading 
information,  internal  decomposition  kinetic  data,  output  interval  specifications 
and  general  program  constants,  nodal  data,  material  property  tables,  pyrolysis 
gas  enthalpy  table,  and  time  dependent  boundary  conditions  table. 

All  this  output  is  fully  labeled  and  is  printed  exactly  as  input  by 
the  user.  Examples  of  these  outputs  are  included  in  Section  4. 


3.2.2 


Surface  Thermochemistry  Tables 


3. 2. 2.1  Edge  Enthalpy  Table 


If  there  is  an  edge  enthalpy  table,  it  is  output  exactly  aS  input. 


3. 2. 2. 2  Surface  Thermochemistry  Table 

This  table  is  output  re-ordered  with  increasing  char  rates,  in  each  , 
section  and  with  a  new  computed  term  in  place  of  the  input  enthalpy  terms.  For 
each  entry  in  the  surface  equilibrium  tables  the  prograin  computes  the'quantity* 


2tw>hi 


I'h  +  mchc  +  mghg  _  qchem  _ 
w  peueCM  peueCM  peueCM 


and  then  outputs  it  as -the  dependent  variable  of  interest  in  the  output  surface 
equilibrium  table.  ^This  term  occurs  directly  in  the  surface  energy  balance, 
Equation- (9)  of  Section  2.4.3  above.  In  the  output,  this  quantity  is  labeled 
CHEM  PROD.  It  has  the  units  (Btu/ftz-sec)/ (lb/ft2-sec)  =  Btu/lb  and  can  be 
thought  of  as  the  chemical  energy  release  per  pound  of  Pe^eCM,  where  pgueCM 
can  be  thought  of  as  the  "scrubbing  flux"  or  "Reynolds  tJux." 

f 

If  edge  tables  are  omitted,  only  the  last  three  terms  are  included. 


!> 


In  the  output  table,  the  quantity  PeuecM  is  abbreviated  to  CM. 

3.2.3  Regular  Output 

3 . 2 . 3 . 1  Introduction 

At  each  output  interval,  as  specified  by  the  user,  the  program  prints 
out  the  current  values  of  ablation  rates,  nodal  temperatures  and  densities,  and 
other  supplementary  information,  as  described  below.  Samples  of  the  regular  out¬ 
put  are  given  in  Section  4 . 


3. 2. 3.2  General  Information 

The  first  line  of  output  shows  the  current  values  of  a  number  of  mis¬ 
cellaneous  quantities: 


Heading 
TIME  STEP 

SURF  ITER 

PROB  OPTN 
SURFACE  RAD  (IN) 

H  WALL 


H  EDGE 

HEAT  COEFF 


CH/CHO 


& 

Total  number  of  computational  cycles  (steps)  required 
to  reach  the  indicated  value  of  output  time.  Note 
that  the  initial  time  is  numbered  as  1. 

Number  of  iterations  required  for  the  surface  energy 
balance  during  the  previous  cycle. 

Current  problem  option  (1,  2,  or  3) 

Current  value  of  surface  radius  for  axisymmetric  geome¬ 
tries;  current  location  of  surface  for  planar  geometries, 
measured  from  original  location. 

Option  1:  Enthalpy  of  frozen  edge  gases  at  the  wall 
temperature . 

Option  2:  Blank 

Option  3:  Enthalpy  of  pyrolysis  gases  at  the  wall  tem¬ 
perature 

The  input  edge  enthalpy  as  determined  by  linear  inter¬ 
polation  in  the  functions-of-time  table. 

The  current  value  of  the  convective  heat  transfer  coef¬ 
ficient  peUeCH  as  linearly  interpolated  in  the  functions- 
of-time  table  and  corrected  for  the  effect  of  transpira¬ 
tion  (blowing)  according  to 

CH/CH  =  C/(e^  -  1)  where  c  =  2Xm/CH  and  CH 
W1  M1  ”l 

is  the  heat  transfer  coefficient  before  being  corrected 
for  blowing;  as  well  as  for  any  radius  ratio  effect. 

The  ratio  C„/C„  ,  indicates  the  amount  of  blowing  cor- 
rection.  1 


3-19 


I 

I 


( 


3. 2. 3. 3 

Ablation  Rates 

The  output 

labels  are  interpreted  as  follows: 

» 

Heading 

Data 

B  PRIME 

The  parameter  B'  =  ®/PeueCM  > 

B  PRIME  G 

\ 

Parameter  B'  based  on  gas  flow  alone:  ma/peueCM 

M  DOT  CHAR 

Current  char  removal  rate  mc,  lb/ft2 sec 

M  DOT  GAS 

Current  pyrolysis  gas  flow  rate  at  the  surface  m  , 
lb/ft2sec.  .  q 

M  CHAR 

Total  char  removal  up  to  current  time,  based  on  a  unit 
area  of  the  original  surface,  lb/ft2. 

M  GAS 

\  Total  pyrolysis  gas  generation  up  to  current  time, 
based  on  a  unit  area  of  the  original  surface,  lb/ft2. 

3.2. 3.4 

Recessions 

and  Recession  Rates  '  ‘  ' 

A  single  output  line  gives  the  current  locations,  as  measured  from 
the  original  location  of  the  surface,  of  the  ablating  surface,  the  char  line 
[defined  as  the  line  of  density  p  +  r  (p_  -  p„)],  and  the  pyrolysis  line,  [de- 
fined  as  the  line  of  density  p  +  r  (p  -  p  )],  as  well  as  the  current  rates  of 

C  P  P  C  M* 

movement  of  these  lines.  The  values  rc  and  rp  are  also  printed  for  convenience. 

3.2. 3.5  Surface  Energy  Flux  Terms' 

Two  lines  of  output  give  the  current  values  of  the  surface  energy  . 

flux  terms  in  Equation  (9)  of  Section  2.4.3,  based  on  one  square  foot  of  the 

present  surface,  and  the  integrated  values  of  these  terms  based  on  one  square 

foot  of  the  original  surface. 

.  (  > 

The  headings  ar.d  data  are  related  as  follows: 


CONVECTED  IN*  ’  qsen  =  PeWHr  -  h^J* 

RADIATED  IN  qra<J  =  ^q^ 


in 

RADIATED  OUT  qrad  =  acT^ 

out 

,  CHEMICAL  .GENERATION  3chem 

.  -CONDUCTION  AWAY  qcond 

*  If  edge  tables  are  omitted,  h€  is  set  to  zero  for  output  purposes. 
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3. 2. 3. 6 


Interior  Energy  Terms 


Two  lines  of  output  show  terms  which  describe  how  the  input  energy" 
of  qcond  is  "accommodated"  or  "partitioned"  in  the  solid  material.  Part  of  it 
is  consumed  in  decomposing  the  plastic,  part  is  consumed  in  sensible  enthalpy 
changes  of  the  solid,  and  part  is  "picked-up"  by  the  pyrolysis  gases  as  they 
pass  through  the  char.  i 

These  terms  are  given  by  the  energy  equation  (3)  of  Section.  2.1.2 
when  these  terms  are  integrated  over  the  total  extent  of  the  ablating  material. 

•  '  The  headings  and  the  terms  are  related  as  follows. 


The  pyrolysis  gas  pick-up  represents  the  energy  picked  up  by  all  the' 
pyrolysis  gas  during  its  passage  through  the  char  layer.  The "decomposition 
absorption  is  simply  the  energy  absorbed  in  isothermally  decomposing  the  virgin 
plastic.  The  loss  at  the  rear  face  is  the  energy  passing  out  the  back  wall  of 
the  ablation  material  (not  the  back  up  materials).  Two  terms  give  the < total 

•  •  a 

rate  of  storage  in  solid  materials.  This  quantity  could  in  fact  be  represented 
by  only  one  term  -- 


A  dx 
A 


s 


‘  i 


but  it  is  more  useful  to  divide  into  the  two^parts  as  indicated.  For  problems 
^with  a  steady  state  temperature  profile  (in  x) ,  the  "Storage  in  Solid"  term 


a 


goes  to  zero.  In  general,  the  relative  sizes  of  the  "Storage  in  Solid"  term 
and  .the  "Convection  with  Solids"  term  indicate  the  nearness  of  the  approach  to 
a  steady  state. 


Heading 

l 

PYROL  GAS  PICK  UP' 


DECOMP  ABSORPTION 


CONVECTION  WITH 
SOLIDS 


Data 


/" 9  SS2),  '  %Vw  +  /  hg  (lf)y  -af 

■/k  '  «  (If) 

J  \  /V 


A  dx 
'Y  As 


-/ipcp  £), 


A  dx 


o  A 
0  s 


3-21 


STORAGE  IN 
SOLID 


LOSS  AT  REAR 
FACE 


H  HI 


The  second  line  of  output  gives  time-integrated  total  values  of  these 
energy  rate  terms,  based  on  one  square  foot  of  the  original  surface. 


3. 2. 3. 7 


Nodal  Data 


/  This  output  gives  the  current  temperature,  density  and  enthalpy*  of 

x~all  the  nodes,  as  well  as  material  number  specifications,  where  material  number 
1  designates  main  material  virgin  plastic,  defined  as  having  a  density  equal  to 
or  greater  than  (pp  -  0.1  lb/fts),  material  number  2  designates  main  material 
i  chari  defined  as* having  a  density  equal  to  or  less  than  (pc  +  0.01  lb/ft3),  ma- 
\  terial  number  0  designates  reacting  char-plastic  mixture  with  intermediate  den- 
! sity,  numbers  21-31  designate  charring  back-up  material  numbers  according  to  the 
table  of  Secjtion  3.1.6,  and  n^terial  numbers  3-10  designate  non-charging  back-up 
materials. 


3. 2. 3. 3  Char  Swell  Information 

If  the  us§r  has  called  for  char  swell  corrections  (Sections  2. 4. 3. 4 

and  3.1.10.2),  two  lines  indicate  the  net  surface  recession  and  current  surface 

'J 

recession  rate  as  affected  by  the  char  swell  effect. 


3.2.4 


Thermocouple  and  Isotherm  Output 


3. 2. 4.1  Listed  Output 

If  thermocouple  and/or  isotherm  output  have  been  called  for,  this 
data. will  be  listed  at  the  end  of  the  regular  output,  giving  thermocouple  tem¬ 
peratures  in  degrees  Rankine  and  isotherm  locations  in  inches  from  the  location 
of  the  original  surface. 

In  both  cases  the  output  data  are  determined  by  a  curve  fitting  rou¬ 
tine  which  extrapolates  beyond  the  physical  boundaries  of  the  system,  so  that 
the  user  must  himself  detect  when  thermocouples  have  ablated  away  and  lost  their 
meaning  and  when  isotherms  are  fictitiously  found  outside  the  material  boundaries. 


Nodal  thermal  conductivity  may  be  output  instead  of  enthalpy.  The  relevant 
flag. punch  for  conductivity  is  in  column  60  of  the  "back  wall  conditions  card" 
described  in  Section  3.1.6  above. 
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3. 2. 4. 2 


Punched  Output 


If  punched  output  of  thermocouple  and  isotherm  data  is  called  for  as 
described  in  Section  3.1.3  above,  the  data  will  be  punched  on  cards  in  the  order: 
time,  surface;  temperature,  thermocouple  x emperatures  (in  the  input  order),  iso¬ 
therm  depths  (in  the  input  order} .  The  format  is  (8F10.4/10X7Fi0.4) . 


3.3.1  . 


DUMPS 

Introductiort 


To  prevent  the  execution  of  computations  that  are  wastefgl  or  probably 
erroneous,  the  program  provides  two  emergency  stops  with  dumps  of  diagnostic 
\  information. 


3.3.2 


Too-Small-Time-Stej 


One  such  dump  occurs  if  the  program  selects  a  time  step  as  small  as 
10  *  seconds.  Computation  ceases,  and  the  program  prints  out  in  a  single  line 
values  of  the  current  time  (which  will,  of  course,  not  in  general  correspond 
to  a  regular  output  time),  the  time  increment  A0,  the  previous  time  increment 
for  comparison,  the  amount  of  problem  time  from  the  previous  time-table  double 
entry  (shift  of  option)  or  discontinuity,  change  in  surface  temperature  during 
last  ‘computation  step,  the  critical  recession  distance  parameter  (equal  to  one- 
fifth  of  the  thickness  of  the  first  node  or  pne-tenthfof  the  minimum  allowable 
last  node  thickness, 'whichever  is  smaller  (A0 , is ' restricted  so  that  recession 
during  the  computation  step  will  not  exceed  this  amount) ,  and  the  surface  recs- 
sion  rate  (ft/sec)  computed  most  recently.  .  f  \ 

Once-  this  data  has  been  output,  the  program  sets  the  problem  time 
exjual  to  the  final  time,  thus  forcing  an  output  of  all  the  current  values  of 
the  standard  output  quantities,  followed  by  a  termination  of  computation. 


3.3.3 


Too-Many-Surface-Iterations  Dump 


The  program  allows  51  iterations  to  find  an  acceptable  surface  energy 
balance.  If  a  balance  has  not  been  obtained  after  51  iteratigns,  the  program 
writes  the  diagnostic  message  "Iteration* Stop"  and  a  block  of  diagnostic  data. 
Since  stops  of  this  kind  are  complex  in  nature,  and  generally  involve  the  sur¬ 
face  thermochemistry  tables  data  deck,  the  dump  data  must  usually  be  communicated 
,  to  the  program  authors  for  analysis. 


{ 


3.3.4  Unacceptable-Surface-Therroochemistry-Table  Stop 

A  series  of  checks  built  into  the  input  routine  serve  to\ detect  com¬ 
mon  errors  in  the  make-up  of  the  surface  thermochemistry  tables.  Discovery  of 
an  error  stops  the  reading  process  and  the  program  prints  out  a  single  line: 

Bad  Surface  Equilibrium  Table  of  Type  - 

Five  error  types  are  detected  as  follows: 


Type 

0  Current  section  of  independent  temperature  entries 

(no  ablation)  is  not  in  descending  order  in  tempera¬ 
ture 

1  Char  rate  in  no-ablation  subsection  has  been  set 
greater  than  a  char  rate  in  the  ablation  subsection. 
Although  the  "zero-char-rate"  can  differ  from  zero, 
they  must  not  exceed  an  actual  ablation  rate  entry. 

2  Edge  table  has  been  omitted  even  though  the  diffusion 
coefficients  are  not  equal  (unequal  diffusion  expo¬ 
nent  /  0) . 

3  Edge  table  has  been  omitted  even  though  C^/C^  /  1. 

4  Inconsistent  unequal  diffusion  exponent  (this  must 
be  uniform  for  all  tables). 

5  Have  called  for  radius  ratio  and/or  burni:  rate  expo¬ 
nent  correction  to  Cjj  although  this  is  i  'an  internal 
radius  geometry.  ° 

3.3.5  cense  Switch  Diagnostic  Pumps 

Sense  switch  calls  built  into  the  Charring  Material  Ablation  Code 
provide  diagnostic  dumps  of  certain  data  during  program  execution.  Usually 
such  dumps  are  of  interest  only  to  the  program  authors,  although  for  reference 
these  dumps  are  summarized  here: 

Sense  Switches  Dump  when  On 


1 


Writes  another  version  of  each  surface 
thermochemistry  table  immediately  after 
the  regular  output  of  this  table;  presents 
final  forms  of  energy  terms  as  actually 
used  in  surface  energy  balance  computations; 
useful  for  analysis  of  sufrace  energy 
balance  failures 

Dumps  nodal  numbers  of  decomposing  back¬ 
ups  ,  after  regular  output  of  input  data, 
for  checking  purposes;  also  dumps  at  each 
standard  output  time  various  density  anu 
pyrolysis  rates  for  nodes  in  charring 
back-up  materials 
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4 


Supplies  standard  output  for  each  time 
step,  regardless  of  formal  output  time 
specifications 

5  Outputs  diagnostic  heat  conduction  quanti¬ 

ties  in  charring  back-ups  each  time  step, 
only  used  for  detailed  analysis 


3 . 4  MISCELLANEOUS 

3,4.1  Running  Time 

Computation  time  depends,  of  course,  on  the  problem  being  computed, 
but  typical  computations  for  charring  materials  run  roughly  in  real  time  on 
the  IBM  7090/7094.  Appendix  B  presents  a  method  for  estimating  run  times. 


3.4.2  Tape  Requirements 

Internal  program  assignments  are  as  follows: 

Tape  Unit  3:  Scratch  used  during  execution 
Tape  Unit. 5 :  Input 
Tape  Unit  6:  Output 
Tape  Unit  7:  Punch 


3.4. 3  Fortran  Deck  Make-Up 

The  Charring  Material  Ablation  program  source  deck  is  in  Fortran  IV 
(or  Fortran  63) . 

The  program  consists  of  the  following  units: 

1 .  Main  Program 

Drives  INP0UT  and  CBM  routines 

2 .  Subroutine  CBM 

Computes  all  time  dependent  results 

3 .  Subroutine*  INP0UT 

Reads  in  all  input  data  and  outputs  same  data  in  front  of 
standard  computed  output 

4 .  Subroutine  LC0UNT 

Counts  lines ,  turns  and  numbers  pages 

5 .  Subroutine  L00K 

Table  look-up  with  linear  interpolation 

6.  Subroutine  SL0PQ 

Quadratic  curve  fit  and  slope  finder 

7 .  Subroutine  0GLE 

Table  look-up  with  cubic  curve  fit 
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8. 


yam».4dya  *■■  u  r  > 


Subroutine  0RDERD 

Ordering  routine 

9 .  Subroutine  0RDERI 

Integer  Array  Ordering  Routine 

10  -  Subroutine  SEQUA 

Orders  according  to  results  of  0RDERD 

11.  Subroutine  SEQUA5 

Orders  up  to  5  Dependent  Variable  Arrays  according  to  results 
of  ORDERD 

3.4.4  Storage  Requirements 

With  "standard"  dimensions  on  subscripted  variables,  the  Charring 
Material  Ablation  program  generally  fits  within  the  32,000  words  of  CDC  and 
Univac  equipment.  The  program  will  not  fit  the  32,000  words  of  the  IBM  7090/ 
7094  machine  by  several  thousand  words.  Simple  dimension  cutting  is  usually 
adequate  to  effect  a  fit  in  most  cases.  In  others,  a  special  overlay  version 
of  the  code  exists,  although  only  minor  improvements  can  be  realized  in  this 
way  since  most  storage  requirements  are  for  data. 

The  first  level  overlay  structure  has  the  following  form: 


Resident  Group 


CMA/CBM  MASTER 

LC0UNT 

L00K 

0GLE 

SL0PQ 


Overlay  #1  Overlay  #2 


INP0UT 
ORDERD 
ORDERI 
SEQUA 
SEQUA 5 


CBM 

THERMS 


3.4.5  Stacking  of  Jobs 

A  series  of  data  decxs  representing  different  problems  may  be  stacked 
for  sequential  computation. 
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SECTION  4 


SAMPLE  PROBLEMS  ' 

This  section  presents  input  listings  and  selected  output  page  list¬ 
ings  for  four  sample  problems  illustrating  various  features  of  the  CMA  Version  3 
code.  The  following  table  summarizes  the  essential  features  of  each  problem. 

As  it  hapoens,  these  examples  are  all  from  rocket  nozzle  technology,  although 
the  CMA  code  is  equally  useful  for  ablating  exterior  heat  shield  analysis,  vari¬ 
ous  combustion  problems,  and,  of  course,  simple  transient  heat  conduction 
studies . 


CMA/CBM  SAMPLE  PROBLEMS 


Sample 

Problem 

Description 

Comment 

No.  1 

Transient  thermal  response  at 
the  throat  location  of  a  carbon 
phenolic  nozzle  eiqpoaed  to  an 
aluminized  solid  propellant. 
Nodal  network  includes  steel 
shell  back-up.  Heat  transfer 
coefficient  is  corrected  for 
radius  change  and  burning  rate 
as  a  function  of  time. 

This  problem  shows  the  effect  of 
burning  rate  correction  on  the 
convective  heat  transfer  coeffi¬ 
cient.  In  addition,  the  input 
of  f-functions  is  demonstrated 
although  only  a  linear  variation 
is  used.  Finally,  since  no  de¬ 
composing  back  up  materials  are 
present,  this  problem  shows  that 
the  CMA/CBM  code  will  do  all  prob¬ 
lems  that  the  CMA,  Version  2  code 
could  do. 

No.  2 

Identical  surface  boundary 
conditions  as  Sample  Problem 

No.  1,  but  the  nodal  network 
contains  seven  silica  phen¬ 
olic  nodes  which  decompose 
when  sufficiently  heated. 

This  problem  demonstrates  the  in¬ 
put  and  output  for  the  CMA/CBM 
solution  of  a  nodal  system  contain¬ 
ing  a  decomposing  back-up  material. 
Also  the  temperatures  of  three 
thermocouples  and  the  depths  of 
two  .isotherms  are  computed  and 
output,  and  the  nodal  therm. a  1  con¬ 
ductivity  is  output  instead  of 
enthalpy. 

No.  3 

Identical  surface  boundary 
conditions  and  nodal  net¬ 
work  as  Sample  Problem  No.  2, 
but  the  radius  in  this  prob¬ 
lem  is  corrected  for  char 
swell. 

The  effect  of  char  swell  on  the 
predicted  recession  is  demon¬ 
strated  in  this  problem. 

No.  4 

Identical  to  Sample  Problem 

No.  2,  except  that  the 
"fissure"  model  is  utilized 
to  describe  pyrolysis  gas 
terms.  (This  requires  a 
change  in  the  surface  thermo¬ 
chemistry  model  and  the  py¬ 
rolysis  gas  enthalpy  table.) 

The  effect  of  "fissure"  model  on 
the  prediction  is  demonstrated. 
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OUTPUT.  SAMPLE  PROBLEM  1  (continued) 
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OUTPUT.  SAMPLE  PROBLEM  1  (continued) 
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OUTPUT.  SAMPLE  PROBLEM  1  (concluded) 


CMA/CBH  ANALYSIS  AT  ABLATIVE  THROAT  Of  A  SOLID  PROPELLANT  POCKET  HOTOP 
TMPOAT  SURFACE  IS  CARBON  PHENOLIC  AT  A  61  OCOPCC  LAYUP  ANGLE*  P-29  ATN 


SILICA  PHENOLIC  DECOMPOSING  BACK-UP  MATL  WITH  f “FUNCTIONS  SAMPLE  NO*  2 


A 

20.25 

0.0  1 

•400*0*  3.00  t .5*0*04 

1000.0 

5 

60.75 

32.40  * 

*480*09  3.00  3.660*04 

600. 

C 

9*. 50 

94*50 

0.0  0.0  0.0 

90000.0  1 

1 91  3  2 

0.0 

120.0 

0.5  2.0  5.0 

1*0  0.3750 

2.0 

10.0 

0.04 

-376.5  0.0  0.0 

-0.345  536. 

1.0  1 

.06 

1.4 

53 j .  1000. 

A 

20.250 

0.0  0 

.140*05  3.0  0.154*05 

tooe. 

8 

60.750 

*0.50  0 

•4*6*10  3.0  0.366*05 

600.0 

C 

129.00 

129.0 

0.0  0.0  0.0 

90000. 

-*407.0  -5260.0  0.0 

-0.315  536.0 

I 

530.0 

0.020 

1.15 

1 

530.0 

0*040 

1 

530.0 

0.040 

l 

530.0 

0.040 

1 

530.* 

0.060 

1 

530.0 

0.100 

1 

530.0 

0.100 

1 

530.0 

0.100 

1 

530.0 

0.150 

1 

530.0 

0.150 

22 

530.0 

0.020 

22 

530.0 

0.0*0 

22 

530.0 

0.040 

22 

530.0 

0.100 

22 

530.0 

0.100 

22 

530.0 

0.100 

22 

530.0 

0.200 

3 

530.0 

0.250 

3 

530.0 

0.250 

0.0 

0.0 

530.0 

1 

530.0 

0.210 

! .040-0* 

0.630  NX4926  .'VlRGfN«60  OCG  LAYUP) 

600.0 

0.360 

2.130-04 

0.630  MX4926  <VI*G1M*60  DEG  LAYUP) 

1160.0 

0.360 

2*470-0* 

0.630  NX4926  (VIRGIN. *0  OCG  LAYUP) 

1500.0 

0.*72 

2.470-04 

0.630  4X4926  (VIRGIN. 6*  OCG  LAYUP) 

2000.0 

*.464 

2.470-0* 

0.630  NX4926  (VIRGIN. 60  OCG  LAYUP) 

3000.0 

#.*«3 

2.470-0* 

0.639  MX4926  (VIRGIN. 60  OCG  LAYUP) 

*000.0 

V.494 

2.470-0* 

0.630  MX4976  (VIRGIN. 40  OCG  LAYUP) 

5000.0 

0.500 

2.470-0* 

0.630  MX4926  (VIRGIN. *•  OCG  LAYUP) 

-1 

6000.0 

0.500 

2.470-0* 

0.630  K.X4926  (VIRGIN. 6*  OCG  LAYUP) 

530.0 

0.210 

2.470-0* 

0.630  NX4926  (CHAR 

•60  OCG  LAYUP) 

1000.0 

0.*30 

7.530-0* 

0.630  NX4926  (CHAR 

•60  OCG  LAYUP) 

1500.0 

0.*72 

2.540-0* 

0.634  MX4926  (CHAR 

•60  OCG  LAYUP) 

2000.0 

*.*«* 

3.240-0* 

6.630  NX4926  (CHAR 

•60  OCG  LAYUP) 

3000.0 

0.*93 

7.190-0* 

0.630  MX4926  (CHAR 

•60  OCG  LAYUP) 

*000.0 

0.*9S 

1.3*0-03 

0.630  NX4926  (CHAR 

•6*  OCG  LAYUP) 

5000.0 

0.500 

2.000-03 

0.634  MX4926  (CHAR 

•60  OCG  LAYUP) 

*1 

6000. 

0.500 

2.750*93 

4.630  NX4926  (CHAR 

*60  OCG  LAYUP) 

530.0 

0.260 

0*0001070 

600.0 

0.275 

0*009109 

1160.0 

0.310 

4.000110 

1500.0 

0.*72 

0.020110 

2000.0 

0**6* 

0.000110 

3000.0 

0.*93 

0.000110 

*000.0 

0.*9J 

O.OOOilO 

-1 

5000.0 

0.500 

0.0001100 

530.0 

0.210 

0.000314 

1000. a 

0.*30 

0.400334 

1500.0 

0.472 

0.000335 

2000.0 

0.*6* 

0.094356 

3000.0 

0.493 

4*440343 

3500.0 

0.495 

0*004395 

*000.0 

0.490 

0*000427 

-1 

5000.0 

0.500 

0*100542 

3 

*•6.00 

500.0 

0.110 

0*00694 

-1 

5000.0 

0.110 

0*00694 

* 

0.0 

0.0 

*•0 

1*0 

-1 

1.0 

1.0 

0.0 

5 

0.001 

0.0 

0.0 

1*0 

0.35 

1.0 

0.600 

0.* 

0.0 

0.90 

0.4 

0.0 

•1 

1.00 

1  0 

0.0 

t 

900.0 

1100.0 

2700.0 

3600.0  4540.0  5*00.0 

63*0.0 

-1702.0 

-9)0.2 

195.6 

2249.0  3*00.0  4696.0 

6152.0 

0.0 

t9*7.0 

904.0 

0.6544  29.2 

10.0 

19*7.0 

904.0 

6.7236  29.2 

HO. *00 

19*7.0 

904.0 

0.7236  29.2 

60.000 

.05 

0. 

1 

120. 

.05 

0. 

0.602 

1.0 

0.625 

29.20 

3500.0 

6.0  914.146  914.106-1 

29.20 

3400.0 

0.0  671.9*5  671.9*5-1 

29.20 

2500.0 

*•*  434.119  434#I|9*| 

29.20 

2000.0 

*.0  2*2.664  202.664-1 

29.20 

1500.0 

4.0  -20.763  -20.763-1 

29.20 

1000.0 

0*0  -231.356  -231.356-1 

29.20 

500.0 

0.4  -405.022  -405.022-1 

29.20*00  0.500*0  •.35*001*30. 72130.  205*. 1*3  2*54.193  1  C*  0. 

29.20*06  6.50606  G. 200003299.77*10.  1697.900  1697.966  l  C*  *. 


INPUT.  SAMPlf  PAOBIOI  2 
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V 


>4.**poo  C.<803fl  0. ISOO* V J9«?4709. 

*.*0050  o.inooftm5.?i06''. 
>9. >0000  0.4*000  0.0^0001050.68700. 

>C-.-">OCO  O.SOOOO  0.  )0O|O?916. 28580. 
?«..>'.<)0*  9.50160  j.ncot<>?soo. ooooo, 
?^,PO.OOO  >,50*00  O.0601O20CO. 00000. 
3&.>rf)0'i  1.60000  O.OO0JO1SOO. 90000. 
>Q.»'-<*9*  O.50OC0  O.OrtOlotooO.OOOOO# 
?9.2,i0O*  0.60000  0,00010  SOO, ooooo, 
•>4.?00OO  0.20000  0.15000 1419.01656. 

*°.?o0oo  o.?oooo  o.p^oeoip^i.iosBo. 

?9.>6000  0.20000  0,1500011 36.66150. 
>9.>«000  0.20000  0. 1*0002981 .10200, 
20.21000  0.20000  0.650002714,33500, 
.'9.2M00  O.ROOOO  O.O10OO26 25, 79420. 

20.26000  0.20000  0.0100071 16.21666, 
>9. >00 oo  0,?0000  0.007501952.70850, 

29.20000  0.20000  9.00500128M.16790* 
?9.?*O0 o  0.20000  0.0001012S2.00290, 
29.2*000  0,20000  o.cooioiooo. 00000. 
20.2600*  0.20000  0.00010  SOO. 00000, 

29.20000  0.10000  0.150001443.87080. 

29.20000  0.10000  0. 200003209. «53140, 
29. >0000  O.JOOOO  3.150001067 .386 10. 

29.20000  G. 1*000  0,100002830.58750* 

29.20000  0.10000  0.0500021  70.40020* 
>9.20000  0.10000  0.040001272.99140, 
29.?0000  9. '0000  0.0 10. 101214.70 890, 

29.20000  0.10000  0.010001131*71570* 

29.20000  0.10000  0.000101095.31176* 

29.20000  0.10000  O.OOOIOIOOO. 00000* 

29.20000  0.10000  O.OOOlC  500. 00000, 
29.20009  0.06000  0.350003446*34560. 

29.20000  0.06000  0.2PC003191 .95550. 

29.20000  0.06000  0. 150001028.23220. 

29.28000  0.06000  0.100002731*50620. 

29.28000  0.06000  0.080002506. 3548C. 

29.20000  0.06000  0 .070002314.54540 . 

29.20000  0.06000  0.065002151.65580. 

29.20000  0.06000  0.040001325.71120. 

79.20000  0.06000  0.050001241.34840. 

29.20000  0.06000  0.030001 151.88190, 

79.28000  0.06000  0*000101018.34140. 

79.20000  0.06000  O.OOOIOIOOO. 00000. 
29.20600  0.06000  0.00010  500*00000. 

20.28000  0.03500  0.350003448.06240. 
29.20006  0.03500  0.2000031 79.22710. 

29.20000  0.03500  0. 1500*7998.60670. 

29.26000  0.03500  0.180002646.98070. 

29.20000  0.03500  0.06000?326.967SO. 

79.20000  0.03500  0.075007166.44010. 

74.20000  0.03500  0.07000J 3J7. 1 7750. 

29.26000  0.01500  0*040001245*53100* 

29.26000  0.03500  0.05606f  1  94.91410. 

29.26000  0.03500  0.00010  899.26150. 

29.20000  0.03500  0.00010  500. 00000. 

29.20000  0.02000  0.350001440*20770. 

29.26000  0.02000  0.2000031 70.86300. 
29*26000  0.02000  0. 150002978. 39280* 

29.20000  0.02000  0* 1 00002582* 76380* 

29.26000  0*02000  0.080002179.65570* 

24.26000  0.02000  0*075001326*44640* 

29.20000  0.02000  0.070001275. 99670. 

29.26000  0.02000  0*06000121 5, Jr 34 60. 

29.20000  0.02000  0.050001170*94960. 
29*26080  0.02000  0,000)0  665.96190* 

29.20000  0.02000  O.OOOlC  600*00000, 

29.20000  0.00010  0.350003450*60590* 

29.26000  0*00010  0.200003156.75570. 

24.20000  0.00010  0. 150002946.10120. 

29.26000  0.00010  0.100002472.02610. 
29*20000  0.00010  0,092502306 *33350* 
29.70000  0.00010  0.08T5021 19.58270. 

29.26000  0.00010  0.085001914.41820. 
?•». 20000  0. 00010  0,082501330.25080. 
?9.2*00*  0.00010  0,080001 299.59080* 
29.76008  0.60010  0 ,06000 1 1 $1 .20 300 . 

29.26000  0.00010  0,050001140.60500. 
79.76006  0.00010  0.00010  600.47750* 
29.76000  0.00010  0.00010  500.00000. 


1557 * 464  1557.984  »  C*  0 

1401-100  1403.100  1  C*  0 

1228.521  1228.521  1  •' *  0 

102«.664  1025.468  1  0  *  0 

485.109  685.109  0  (  0 

298.69?  *98.89?  0  0*  0 

-256.818  -250.818  0  31203*  0 

-911.760  -933.760  0  .1.20  «•  C 

-1305. 31t>-l  305.336  0  AI.20J*  0 

18(8.674  1418.079  1  i  •  0 

1  142.744  1  142.244  1  *  ♦  0 

I  1 42.64  7  1  142.947  l  C*  0 

907. *14  907.434  I  O  0 

602.751  402. 7S3  1  C*  0 

435.126  435.126  1  C*  0 

145.414  145.414  |  C*  0. 

47.720  42.720  1  C*  0 

-494.248  -494,246  l  C*  0 

-552.716  -552.7J8  1  C-  0 

-1012.490-1012.490  0  AL203*  0 

-1353.107-1353.107  0  AL203*  0 

1718.127  1718.127  1  C*  0 

1182*973  1182.973  1  C*  0 

949.446  949.646  1  C*  0 

656.668  656*888  1  C*  0 

138.870  l 36,870  1  C*  0 

-526.689  -526.669  l  C*  0 

-629.431  -629.431  1  <*•  0 

-799.609  -799.609  l  C*  0 

-664.688  -684.668  If*  0 

-1045* 958-1045. 956  0  8:203*  0 

-1374.627-1374.627  0  At  203*  0 

1674*103  1674*303  1  C*'  0 

1111*021  1 t 11*021  1  C*  0 

659.936  659.936  l  C*  0 

531.122  531.322  1  C*  0 

344*691  144.691  1  C*  0 

211.524  211.524  1  C*  0 

109.406  109.606  1  C*  0 

-461.916  -461.916  1  C*  0 

-561.601  -563.801  1  C*  0 

-756.859  -756.659  1  C*  0 

-1033.189-1033.389  1  C*  0 

-1061.054-1061.054  0  41203*  0 

-1364.661-1304.661  0  AL203*  0 

1645.676  1645.676  1  C*  0 

1063.164  1063.164  1  C*  0 

799.242  799.242  I  C*  0 

440.337  440.317  1  C*  0 

206.774  206.774  1  C*  0 

106.992  106.99?  1  C*  0 

-452.619  -452.619  1  C*  0 

-578.870  -578.670  1  C*  0 

-670.121  -670.321  1  C*  0 

-1195.57I-U95.S71  1  C*  0 

-1191-193-1191*193  0  C*  0 

1626.023  1628.023  1  C*  0 

1033.120  1013.120  1  C*  0 

760.609  '60.609  1  C*  0 

378.625  376.825  l  C*  0 

76.164  76.164  |  C*  0 

-465.197  -465. 197  1  C*  0 

-533.272  -533.27?  1  C*  0 

-632.566  -632.586  1  C*  0 

-716.460  -716.460  1  C*  0 

-1350. 252-1 158.252  1  C*  0 

-1395.265-1395.265  0  AL203*  0 

1604.023  1604.023  1  C*  0 

992.289  99?. *80  l  C*  0 

707.221  70T.221  l  C*  0 

284.658  284.656  1  C*  0 

175.096  175.096  1  C*  0 

62.729  62.729  1  C*  0 

-56.221  -56.221  1  C*  0 

-466.518  -466.516  l  C*  0 

-503.670  -503.670  I  C*  0 

*696.006  -698.006  l  C*  0 

-780.959  -780.959  I  C*  0 

-1384.243-1384.243  1  C*  0 

-1400. 855-1400. 855  0  A1201*  0 


IHPUr,,  SAMPLE  PROBLEM  2  (cohcluded) 
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APPENDIX  A 


SOME  SUGGESTIONS  FOR  SELECTING  NODAL  SIZES 

A.  1  GENERAL  REMARKS 

Experience  with  particular  problems  will  always  be  the  best  guide  in 
selecting  the  sizes  of  the  finite  difference  nodes  within  the  heated  material. 
Nevertheless,  it  is  possible  to  set  down  some  very  rough  guidelines  for  nodal 
size  which  require  only  a  preliminary  estimate  of  the  overall  character  of  the 
solution  history  rat)*"  *han  detailed  experience  with  closely  related  problems. 
The  following  subsection  describe  how  such  preliminary  estimates  cart  be  turned 
into  moderately  useful  nodal  size  assessments. 

A. 2  NODAL  SIZES  RELATED  TO  TIME  CONSTANT  INFORMATION 

Most  transient  problems  have  some  identifiable  time  constant  which 
can  be  used  to  select  a  nodal  size.  This  time  constant  is  either  the  basic  re¬ 
sistance  times  capacitance  time  constant  of  the  system  (usually,  of  course,  the 
chief  resistance  is  found  in  the  ablating  material  itself,  since  the  outer  heat 
transfer  conductance  is  very  high  in  most  ablation  problems)  or  some  time  con¬ 
stant  associated  with  a  transient  boundary  condition.  The  second  case  is  the 
simpler  of  the  two  to  visualize:  obviously,  we  wa\it  the  time  constant  of  the 
nodal  zones  to  not  be  significantly  larger  than  the  time  constant  of  the  bound¬ 
ary  condition,  0cb.  The  time  constant  of  a  node  is  simply 
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where  p,  Cp,  and  k  are  evaluated  at  the  nodal  temperature.  Thus,  the  first 
nodal  size  criterion  is  simply 
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Fox  problems  with  relatively  steady  boundary  conditions,  Equation 
(A-3)  will  not  yield  useful  results  and  one  must  consider  the  "natural"  time 
constant  of  the  ablating  material,  i.e.,  the  first  case  mentioned  above.  For 
a  problem  with  negligible  recession  and  steady  boundary  conditions,  no  single 
set  of  nodes  can  be  appropriate  for  the  entire  problem  history.  If  we  let  0^ 
stand  for  the  time  around  which  interest  centers,  then  the  standard  similarity 
so3.ution  for  the  transient  temperature  in  a  constant  properties  slab  initially 
at  constant  temperature  T^  with  steady  surface  temperature  Tq  impressed  at 
0  =  0  shows  that  the  dimensionless  temperature  (T  -  T.)/(T  -  T. )  has  decayed 
to  less  than  0.1  at  x/2(a0^)^  =  2.  If  we  want  10  or  so  nodes  distributed  over 
this  region,  then 

x  w 
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which  naturally  enough  is  similar  to  the  result  of  Equation  (A-3) ,  ^ven  though 
two  different  time  concepts  are  involved. 


A. 3  NODAL  SIZES  RELATED  TO  SURFACE  RECESSION  RATE  INFORMATION 

If  (as  is  usual)  the  problem  has  surface  recession  S,  then  some  addi- 
tional  time  step  considerations  must  be  made.  If  S  and  Tw  are  constant  values 
impressed  at  time  zero  for  a  constant  properties  slab  initially  at  T^,  the  simple 
analysis  shows  that  the  in-depth  temperature  profile  as  measured  from  the  heated 
surface  ultimately  attains  a  steady  shape 


(A- 5) 


This  profile  decays  to  0.1  for  Sx/a  *  3;  if  we  want,  say,  10  nodes  in  this  in¬ 
terval  then 

5  =  ]^  =  ~  =  ° .  3«/S 
10S 

v/ 

However,  a  slightly  tighter  criterion  than  this  derives  from  a  basic  accuracy 
study  of  the  particular  finite  difference  equations  used  in  CMA.  This  study, 
reported  in  Appendix  B  of  Reference  A-l  below,  indicated' that  for  good  in-depth 
accuracy  one  should  satisfy  the  restriction 

6  ~  O.la/S 
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A  final  result  of  interest  may  be  obtained  by  noting  that  the  con¬ 
stant  S,  constant  Tw  problem  described  here  does  have  a  transient  beginning 
and  approaches  the  steady  state  only  as  a  limit.  References  A-2  and  A-3  show 
that  this  limit,  as  measured  by  the  closeness  of  approach  of  Qcon(j  to  its 
steady  state  value  is  reached  to  within  a  few  percent  at 

0  *  a/S2  (A- 8) 

It  is  very  interesting  to  observe  that  if  we  consider  the  time  given  by  Equa¬ 
tion  (A-8)  as  the  "time  of  interest"  and  substitute  it  into  the  time  constant 
result  (A-4),  we  obtain  the  result 


6  *  0 . 4a/S  (A-9) 

which  harmonizes  rather  remarkably  with  Equation  (A-6),  derived  from  steady 
state  profile  considerations. 

A. 4  CONSEQUENCES  OF  SELECTING  INAPPROPRIATE  NODE  SIZES 

The  selection  of  node  sizes  much  smaller  than  necessary  generally 
merely  increases  the  computing  time  to  far  above  what  could  be  achieved  with 
more  reasonable  nodal  sizes.  In  theory,  extremely  small  nodal  size  could  intro¬ 
duce  excessive  numerical  round-off  error  problems,  although  this  has  never  been 

* 

observed  in  practice  with  the  CMA  program. 

The  choice  of  too  large  nodes  usually  results  in  a  non-obvious  lack 
of  accuracy  for  the  entire  prediction,  especially  for  the  in-depth  temperature 
profile.  Important  problems  should  be  computed  with  two  different  choices  of 
nodal  sizes  for  a  check  of  the  accuracy  attained.  Gross  violations  of  restric¬ 
tion  (A-7)  sometimes  cause  drastic  and  obvious  oscillations  in  the  in-depth 
predictions . 

A. 5  REMARKS  ON  NODELETS 

As  noted  in  Section  2.1.1  above,  the  pyrolysis  events  in  the  CMA  pro¬ 
gram  are  computed  in  a  grid  of  "nodelets"  finer  than  the  nodal  grid  used  for  the 
energy  calculations.  The  number  of  "nodelets"  per  node  may  be  selected  by  the 
user.  Experience  to  date  has  not  been  sufficient  to  lay  down  any  general  rules 
for  selecting  the  number  of  nodelets  per  node  in  various  applications.  Ten 
appears  sufficient,  although  many  problems  seem  to  do  well  with  fewer.  The  num¬ 
ber  of  nodelets  per  node  often  has  an  appreciable  effect  on  computation  time 
(see  Appendix  B  below) ,  and  for  an  extensive  series  of  calculations  it  is  worth¬ 
while  to  minimize  this  number  (by  experiment)  to  reduce  costs. 


SUMMARY 
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The  user  confronted  with  an  ablation  problem  outside  the  range  of 
his  prior  experience  can  draw  upon  several  rough  guides  to  select  nodal  sizes. 
Locally  throughout  the  material,  nodal  thicknesses  should  satisfy  the  basic 
thermal  profile  penetration  relation 

<S  *  O.4{a0i),s  (A- 10 ) 

where  0^  is  a  representative  real  time  of  principal  interest  and  a  is  evaluated 
locally  throughout  the  material.  For  example,  suppose  the  user  is  studying 
nozzle  ablation  during  a  60  second  rocket  firing  and  specifies  output  every  10 
seconds.  The  "real  time  of  interest"  to  the  user  is  evidently  about  10  seconds, 
since  the  first  output  will  occur  then.  The  user  would  set  0^  equal  to  10  sec¬ 
onds  in  Equation  (A-10) ,  and  select  a  values  for  the  char  for  estimating  nodal 
spacing  near  the  surface.  Values  of  a  for  the  virgin  plastic  would  be  used  for 
estimating  the  nodal  spacing  deeper  in  the  material. 

1 

For  time  varying  boundary  conditions  with  a  "time  constant"  or  rime 

scale  of  0  .  >  i 

cb  x 

6  *  (a6cb)h  (A- 11) 

if  this  variation  is  to  be  "followed"  accurately  by  the  in-depth  solution.  For 
example,  suppose  the  user  is  analyzing  heat  shield  ablation  during  a  reentry  for 
which  the  "heating  pulse"  rate  rises  from  zero  to  a  maximum  in  60  seconds  and 
then  decays  again  in  an  additional  100  seconds.  Speaking  in  very  approximate 
terms,  the  "time  constant"  of  such  a  pulse  represents  the  time  required  for  the 
pulse  to  rise  about  half  way  from  "trough"  to  peak  or  to  sink  roughly  half  way 
from  peak  to  trough.  This  gives  time  constants  of  30  and  50  seconds;  the  user 
would  naturally  select  the  shorter  one  for  use  in  Equation  (A-ll) . 

Finally,  if  recession  is  important,  nodal  sizes  should  also  satisfy 

6  *  O.lo/S  (A-12) 

to  preserve  accuracy  of  the  solution. 

Most  generally  users  prefer  to  grade  the  nodal  sizes  within  the  gen¬ 
eral  limits  set  by  these  guidelines,  beginning  with  smaller  nodes  near  the 
heated  surface  and  expanding  to  rather  large  nodes  (even  violating  these  guide¬ 
lines)  in  regions  where  the  thermal  pulse  does  not  penetrate,  or  just  begins  to 
penetrate.  It  is  recommended  that  the  user  make  the  first  node  half  the  thick¬ 
ness  of  the  second  node  in  all  cases,  since  the  program  assumes  the  first  node 
is  in  fact  a  half  node  with  only  half  the  usual  number  of  nodelets. 
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APPENDIX  B 

ESTIMATION  OF  EXECUTION  TIME 


Estimates  of  computer  run  time  can  generally  be  made  by  determining 
the  required  computation  time  per  time  step  and  multiplying  this  by  the  esti¬ 
mated  number  of  time  steps  needed  (equal  to  problem  time  divided  by  time  step 
size) .  The  following  sections  will  discuss  this  approach.  A  quicker  estimated 
adequate  for  some  purposes,  may  be  obtained  from  the  simple  empirical  statement 
that  most  CMA  problems  consume  computer  time  at  a  rate  of  about  20  to  30  percent 
of  the  real  problem  time.  This  estimate  applies  to  charring  problems  of  about 
50  nodes  with  recession  rates  less  them  about  25  mils/sec. 

B.l  ROUGH  ESTIMATES  OF  COMPUTER  TIME  PER  TIME  STEP 

"Typical"  problems  of  the  type  cited  above  have  been  observed  to 

t 

require  about  50,000  usee  to  100,000  ysec  per  time  step  on  Univac  1108  and  CDC- 
6600  machines,  and  about  500,000  usee* per  time  step  on  the  IBM  7094  machines. 

B.2  '  REFINED  ESTIMATES  OF  COMPUTER  TIME  PER  TIME  STEP 

A. refined  estimate  of  the  computer  time  required  per  time  step  can 
be  obtained  from  a  count  of  arithmetic  operations.  The  following  table  of  the 
number  of  operations  in  each  section  of  the  Charring  Material  Ablation  program 
provides  a  useful  estimate  of  execution  time  as  a  function  of  the  number  of 
nodes,  the  number  of  nodelets  per  node,  and  so  on. 

The  input  subroutine  and  all  other  operations,  such  as  output,  not 
performed  each  time  step,  'have  not  been  included.  Only  floating  arithmetic  is 
considered . 


As  an  example,  with  the  average  operating  times  for  the  IBM  7094 
(add  and  subtract  =  14  s,  multiply  =  7  s,  divide  =  12  s,  expf  =  188  logf  = 
226  s) ,  this  table  yields  the  following  floating  arithmetic  time  esem  .te  per 
time  step. 


TABLE  OP  FLOATING  OPERATIONS 
FOR  ONE  TIME  STEP 


Add  . 

Subtract 

Multiply 

Divide 

Expf 

Logf 

.  Miscellaneous 

Nodal  Computa¬ 
tions,  per  Node 

32 

36 

64 

25 

0 

0 

1 

General  Decompo¬ 
sition,  per  De¬ 
composing  Node 

22 

15 

13 

3 

0 

0 

Component  Decompo¬ 
sition,  per  Com¬ 
ponent  for  each 
Decomposing 
Nodelet 

1 

3 

4 

3 

3 

2 

Surface  Calcula¬ 
tions,  per 
Iteration 

25 

39 

50 

10 

1 

0 

T  «  [(1136  •  I  +  645) J  *  F  +  1700]N  +  K(1534)  psec/step 


where 

T  *  execution  tine  for  a  single  tine  step  (psec) 

I  *  number  of  decomposing  components  in  the  material. 

J  =  number  of  node lets  per  node 
F  *  fraction  of  nodes  decomposing 
N  =  total  number  of  nodes 

K  =  number  of  iterations  in  surface  energy  balance 

\  • 

In  this  equation,  F  should  be  an  average  value  for  the  fraction  of  nodes  actively 
decomposing  during  the  problem  (usually  between  1/3  and  1/2) .  A  study  of  numer¬ 
ous  problems  reveals  that  K  averages  very  close  to  3, 

*  The  resulting  total  time  must  be  multiplied  by  a  factor  to  account 

for  "administrative"  arithmetic  such  as  DO  loop  indexing.  Experience  with  this 
program  on  the  7094  suggests  a  factor  between  1.4  and  1.5 

B. 3  DETERMINATION  OF  NUMBER  OF  TIME  STEPS  REQUIRED 

'  J 

To  estimate  total  execution  time,  the'  time  per  time  'step  T  may  be 
multiplied  by  the  total  number  of  time  steps,  a  number  approximately  determined 
by  the  user  in  the  choice  of  the  maximum  time  step  allowed  (see  Section  3.1.3 
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above) .  For  some  problems,  however,  a  built-in  limit  based  on  the  recession 
rate^ turns  out  to  be  smaller.  This  limit  prevents  the  time  step  A0  from  allow¬ 
ing  more  than  the  thickness  of  one  nodelet  to  ablate  during  a  time  step;  thus 

A0  =  min  (26^,62,63, <$4, . .  .6NL  }/SJ 

min 

■> 

where  the  6’s  are  the  nodal  thickness  in  the  main  charring  material  and  6NLm.n 
is  the  minimum  allowed  thickness  for  the  last  ablating  node  (Section  3.1.3, 
p.  3  *3  above) . 

B.  4  OTHER  CONSIDERATIONS 

Additional  machine  time  will  be  consumed  by  input  and  loading  opera¬ 
tions,  overlay  or  special  tape  handling  (if  necessary) ,  and  output  operations, 
the  details  of  which  cure  peculiar  to  individual  computing  facilities. 
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